Memoirs of the Faculty of Science, Kyoto University, Series of Physics, Astrophysics, 
Geophysics and Chemistry, VoLXXXXIV, No.l, Article 2, 2003 



J /ij) Production in p+p Collisions 
at ^5 = 200 GeV 

By 

Hiroki Sato 

Department of Physics, Faculty of Science, Kyoto University 
Kyoto 606-8052, Japan 

( Received December 9, 2002) 
Abstract 

Total and differential cross sections for inclusive J/tj: production 
have been measured in ^/s = 200 GeV p+p collisions at Relativistic 
Heavy Ion Collider in Brookhaven National Laboratory. J /ip particles 
have been clearly identified via fJ,'^fJ,~ decays measured in the forward 
muon spectrometer covering 1.2 < y < 2.2 and via e~^e~ decays mea- 
sured in the mid-rapidity spectrometers of the PHENIX experiment. 
Details of the muon channel measurement are presented in this paper 
based upon an 81-nb~^ integrated luminosity obtained in the second 
run period at RHIC. 

Corrected with the detection efficiency, the branching fraction times 
rapidity-differential cross section at forward rapidity Br{J/ilj 
daj/^/dy\y=i,'r = 37 + 7 (stat.) ± 11 (syst.) nb has been obtained. 
The average transverse momentum of the J/'^'s, {pt) = 1-66 it 0.18 
(stat.) lb 0.12 (syst.) GeV/c, is slightly higher than lower-energy re- 
sults and consistent with a model expectation based on perturbative 
QCD. The total production cross section ctj/^ = 3.8 it 0.6 (stat.) it 
1.3 (syst.) /xb has been extracted by fitting both the muon and electron 
channel measurements, which is related with lower-energy results by 
gluon distribution functions currently available. Its normalization is 
well reproduced by both the color-octet model and color-evaporation 
model. 

The results have established the baseline necessary for the QGP 
search in Au+Au collisions as well as gluon polarization measurements 
in polarized p+p collisions at RHIC, where the unpolarized cross sec- 
tions presented here will determine the absolute normalization. 



2 



J/V' Production in p+p Collisions at ^/s = 200 GeV 



Contents 



1 Introduction 4 

1.1 Heavy ion pliysics 4 

1.2 Spin pliysics 5 

1.3 J/ip production in hadron-hadron collisions 7 

2 Experimental setup and data acquisition 13 

2.1 The RHIC accelerator complex 13 

2.2 Coordinates and formulae 14 

2.3 The PHENIX experiment overview 15 

2.4 Interaction trigger counters 16 

2.4.1 The Beam-Beam Counters 17 

2.4.2 The Normalization Trigger Counters 17 

2.5 The Muon Arms 18 

2.5.1 The Muon Tracker 19 

2.5.2 The Muon Identifier 25 

2.5.3 The MulD NlM-logic trigger for LVL-1 30 

2.6 The PHENIX DAQ system 31 

2.7 Global Level- 1 triggers 32 

2.7.1 Minimum bias trigger 32 

2.7.2 Muon triggers 33 

2.8 Onhne Monitoring 34 

2.9 Run and data summary 35 

3 Data analysis 37 

3.1 Data sample 37 

3.2 Muon track reconstruction 37 

3.2.1 MulD road finding 38 

3.2.2 MuTr cluster finding / fitting 38 

3.2.3 Muon track finding / fitting 39 

3.3 Single muons 41 

3.3.1 Event vertex distribution 41 

3.3.2 Track-road matching 43 

3.4 Dimuon mass and Nj/^ 43 

3.5 J/ip detection efficiency 47 

3.5.1 Detector acceptance 49 

3.5.2 MulD efficiencies 50 

3.5.3 MuTr efficiencies 55 

3.5.4 BBC efficiency for p+p — > J/ipX events 59 

3.5.5 Total detection efficiency 59 



H. D. Sato 3 

3.6 Integrated luminosity 60 

3.6.1 Number of minimum bias triggers 61 

3.6.2 p+p inelastic cross section 62 

3.6.3 BBC efficiencies 62 

3.7 Average pt 67 

3.7.1 Uncorrected (pt) 67 

3.7.2 Efficiency corrections 68 

3.7.3 Momentum offset and resolution correction 68 

4 Results and Discussions 69 

4.1 Differential and total cross sections 69 

4.1.1 Transverse momentum distribution 70 

4.1.2 Rapidity distribution and gluon distribution 71 

4.1.3 Total cross section and its comparison with lower-energy 
results 76 

4.2 Average pr 79 

4.3 Future measurements 82 

5 Conclusion 87 
A Estimation of veirious contributions to the inclusive J/ip yield 88 
B Attached figures' captions 90 



4 



J/V' Production in p+p Collisions at ^/s = 200 GeV 



1 Introduction 

Since the sensational discovery of the J/ip particle [1,2] as a triumph 
of the Quark Model, heavy- quarkonium production in hadron-hadron col- 
lisions and in photon-hadron collisions has drawn attentions because it in- 
volves both perturbative and non-perturbative aspects of Quantum Chromo- 
dynamics (QCD), which describes interactions between quarks and gluons. 
A sufficiently large scale Q ~ Mq (heavy quark mass) enables perturbative 
treatments of the production of a heavy-quark pair. On the other hand, 
its hadronization into a quarkonium state, with a scale of a typical hadron 
size (~ Aqcd), is considered to be a non-perturbative phenomenon, where 
perturbative calculation is not applicable. 

Thanks to narrow resonances found in the invariant mass of lepton pairs 
and relatively large cross-sections especially for the case of J/ip (IS), abun- 
dant experimental data have been accumulated on heavy quarkonia. Results 
of cross sections, polarization, relative yields such as the ip' to J/ip ratio, 
have led to better understanding of quarkonium production. However the 
production mechanism involving non-pcrturbativc QCD phenomena is still 
controversial. To understand it, more data on various observables in a wider 
energy range are needed. 

We have added new measurements of J/ ip in p-|-p collisions at unexploited 
energy (y^ = 200 GeV) with the PHENIX experiment at Rclativistic Heavy 
Ion Collider (RHIC). The results are expected to provide better understand- 
ing of non-perturbative aspects of the hadro-production of J/ip, as well as 
further knowledge of perturbative QCD. 

In addition, our measurements play crucial roles for both two innovative 
measurements realized for the first time at RHIC, heavy ion physics and 
spin physics, which will bring new insights into both perturbative and non- 
perturbative aspects of QCD. 

1.1 Heavy ion physics 

Provided nucleons are in high temperature and/or high density environ- 
ment, hadronic matter is expected to turn into the Quark-Gluon Plasma 
(QGP), a new phase of matter whose degrees of freedom are quarks and glu- 
ons instead of hadrons. Understanding the behavior of bulk matter governed 
by QCD elementary degrees of freedom and interactions, and studying how 
it turns into hadronic matter, offers challenging perspectives and touches 
fundamental issues in the study of QCD in its non perturbative regime, such 
as the nature of confinement and chiral symmetry breaking. 

RHIC is the first heavy-ion coUider in the world to create the highest tem- 
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perature environment experimentally. In the most central Au+Au collisions 
at the center-of-mass energy ^/s = 200 GeV per nucleon-nucleon collision, 
initial temperature, Tq, is expected to reach as high as Tq ~ 0.5 GeV [3] at 
which the formation of QGP is expected. 

A heavy quarkonium, a bound state of a heavy quark (Q) and anti heavy 
quark (Q), is detected via a leptonic decay with the PHENIX detector at 
RHIC. It is considered to be one of the best probes for the earliest stages of 
collisions since QQ pairs can be produced only in hard parton interactions 
(dominantly gluon fusion at RHIC energies) and final state leptons do not 
suffer from interactions with hadronic matter in the later stages of the colli- 
sions. In deconfined (QGP) matter, Q and Q tend not to form a quarkonium 
and diffuse away from each other instead, because attractive color forces be- 
tween them are reduced by a Dcbye-type screening. Therefore, suppression 
of quarkonium yield is expected in QGP [4]. Another theoretical model pre- 
dicts enhancement of quarkonium production due to recombination of Q and 
Q produced in different nucleon-nucleon collisions because of their higher 
mobility inside the QGP [5]. If we observed an anomalous suppression or 
enhancement of quarkonium production in Au+Au collisions, it could be in- 
terpreted as evidence of QGP. This study definitely requires the reference of 
unsuppressed cross section, which can be best served by p+p collisions at the 
same energy. There has been, however, no data of J/ip production in p-|-p 
collisions at -y/s = 200 GeV, and such measurements play a crucial role to 
establish the basis of the QGP search in heavy ion collisions. 

The J/ip has the largest production cross-section among all the heavy 
quarkonia and possible to be measured even with such a small integrated 
luminosity as obtained in the RHIC Run-2 period (2001-2002) where acceler- 
ator performance was in the stage of progress. With an increased luminosity, 
measurements and comparison with other charmonium and bottomonium 
states are important, since different binding radii for each quarkonium will 
result in different degrees of suppressions. 

1.2 Spin physics 

Spin is one of the fundamental properties indigenous to nuclei and ele- 
mentary particles. Studies of spin have been providing deep understandings 
of properties of particles and their interactions. 

An example is the studies of nucleon spin structure using deep inelas- 
tic scattering (DIS) of polarized leptons from polarized nucleons. Since the 
1970s, the polarized structure functions of nucleons, gi^'^\x) have been mea- 
sured using polarized electrons or muons and polarized nuclear targets [6]. 
As a result, the net quark polarization inside a proton AE has turned out to 
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be as small as 0.1 to 0.3 which is much smaller than the naive expectation 
(~ 1). This discrepancy is called the "proton spin puzzle". A significant frac- 
tion of the proton spin is allegedly carried by gluons in analogy to the case of 
momentum. However, polarized gluon distribution Ag{x) = g{x)^ — g{x)~ 
where g{x)^{g{x)~) is distribution of gluons polarized in the same (opposite) 
direction as the proton helicity, suffers from large theoretical uncertainties in 
extraction from the DIS results since gluons do not couple to photons at the 
leading order. 

RHIC has been added a new feature of accelerating polarized protons up 
to 250 GeV/c momentum keeping 70% polarization with the advent of inno- 
vative technologies such as Siberian Snakes [7] to keep proton polarization 
through acceleration. Studying polarized phenomena at the perturbative 
QCD regime is now accessible using polarized p-|-p collisions at the highest 
energies up to ^/s = 500 GeV [8], which is much higher than the previous 
polarized p-|-p experiment {y/s = 19 GeV) [9]. For example, we plan the 
first direct measurement of Ag{x), through double-spin asymmetries {An) 
for cross sections of gluon-initiated subprocesses such as the direct-photon 
production [10]. It is important to measure An for various subprocesses 
to elude theoretical uncertainties to extract Ag{x). Quarkonium production 
is one of them. An asymmetry A^^ for quarkonium ipg production can be 
written as 

_ Ag{xi) Ag{x2) ^gg^^gX 

- g(x,) g{x,) 

where xi{2) denotes a momentum fraction of each parton and a^^^''^ is 
a partonic asymmetry for the subprocess g + g ^ ipgX^. An advantage 
of this channel is an experimental feasibility due to a relatively large cross- 
section and small background especially for the case of J/ ip as already stated. 
Also the subprocess (gluon fusion) can be identified with a small (< 10%) 
contamination from others. Therefore we will be able to measure Aj^^ for 
g + g —>■ J/ ipX with a small uncertainty. 

However, the extraction of Ag[x) from A^'^ suffers from large theoret- 
ical ambiguities due to the lack of knowledge of a^^^^, which is sensitive 
to the production mechanism. Therefore the elucidation of the production 
mechanism is a crucial key to extract Ag{x) from Aj^^ . 



In summary, measurements of J/ip in p-|-p collisions at RHIC are vital 
to: (1) discuss the formation of QGP in Au-|-Au collisions and (2) clarify 
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the mechanism of quarkonium production to obtain A.g{x) in polarized p+p 
coUisions. 

1.3 J/ip production in hadron-hadron collisions 

Current theoretical and experimental knowledge of hadro-production of 
charmonia is described in the following. 

Production of J/ijj and other charmonia in hadron-hadron collisions is 
understood in the following framework based upon the QCD factorization 
theorem [11]. Figure 1 shows an example of Fcynman diagram for hadro- 
production of J/ip. According to the factorization theorem, the cross section 
to produce a charmonium -0 in a collision of hadron A and B, a{AB — > ^jJX) 
can be factorized into: 

• fa/A{x,Q): Probability for a parton a to be found in a hadron A, 
called a Parton Distribution Function (PDF), where x is a momentum 
fraction of a to A and Q is a scale for the parton interaction, 

• fb/B{x, Q): Probability for a parton b to be found in a hadron B and 

• a{ab 'ipX): Cross section for the partonic subprocess ab ^X, 
where ip is produced in the hard scattering of the parton a and 6, 

and written as 



a{AB i/jX) = ^ dxidx2fa/A{xi,Q)fb/B{x2,Q)a{ab ipX). 



Following 2— >1 subprocesses contribute to low-pr (transverse momentum) 
production of -0 (< 1 GeV/c): 

where q and g denote a quark and a gluon respectively. Figure 1 is an example 
of the g + g ^ ip subprocess. Medium and high pt ip^s are produced in the 
following 2— >2 subprocesses: 



ab 




• 9 + 9 

• 9 + q 

• q + q 
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Figure 1: An example of J/ip hadro-production. In this example, a J/ip 
meson is produced in the subprocess g + g ^ J/ip in a. p+p collision. The 
oval between the cc pair and the J/ip represents the formation of the J/ip 
from the cc pair which is a non-perturbative phenomenon. 

At RHIC energies (^i = 200 to 500 GeV) g + g and g + q subprocesses 
contribute to cross sections significantly. 

Production of from initial partons is further separated into two steps: 
production of a cc pair in the hard scattering of the initial partons (step 1) 
and hadronization of a charmonium from the cc pair (step 2). This is pos- 
sible since energy scales for each step are well separated, that is, 2Mc ~ 3 GeV 
(Mc is the charm quark mass) for step 1 and Aqcd ~ 0.2 GeV for step 2 
which is the inverse of a typical hadron size. The cross section for step 1 
can be calculated using perturbative QCD. On the other hand, the proba- 
bility for step 2 is not calculable with it, because it is a non-perturbative 
phenomenon. There are some theoretical models for charmonium production 
each of which assumes a different assumption and treatment on step 2. 

The color-evaporation model (GEM) or the semi-local duality approach [12, 
13], born in the 1970s, simply ignores the color and other quantum numbers 
of cc pairs and assumes a certain fraction of them, which is needed to be 
determined from experimental data, forms each charmonium state through 
multiple soft-gluon emissions as illustrated in Fig. 2. The GEM describes 
experimental data on low-p^- or p^-integrated results well, where the soft 
gluon picture is expected to be valid. For example, the GEM describes J/ip^s 
total cross sections in both hadro-production and photo-production at lower 
energies [14]. Also the GEM prediction of zero polarization (spin-alignment) 
of J/ip^s is consistent with the lower-energy experiments where low-p^ con- 
tribution is dominant [43], but contradicts to the GDF data at medium and 
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highpT [47]. 




Figure 2: A schematic diagram for J/ip production from gluon fusion with 
the color-evaporation model. Incoming protons are omitted in the figure. 
Multiple soft-gluon emissions destroy the information on quantum numbers 
of the cc pair. 

A more sophisticated model born in the 1980s, the color-singlet model 
(CSM) [15, 16], requires a cc pair to be the color-singlet state and have the 
same quantum numbers as the charmonium to be formed. Figure 3 shows 
an example of the lowest order production of a J/ i/j with the CSM where 
the cc pair should be in ^'^'^^Lj = ^Si and the color-singlet state as the 
J/ip. It should be noted that an additional hard- gluon emission is necessary 
to conserve the C-parity. The CSM can unambiguously predict production 
cross sections for each charmonium without any free parameters and has 
explained pr distributions of J/ip production at ISR energies {y/s = 30 to 63 
GeV) reasonably well [15]. However the CSM failed to explain pt differential 
cross sections of the Tevatron data in p-|-p collisions at y/s = 1.8 TeV by 
large factors (30 ~ 50) [45]. CSM predictions do not agree either with the 
total cross sections at lower energies by a factor of about 20 [22] . 

To explain these discrepancies, the color-octet model (COM) [17] has been 
developed in the 1990s based upon the non-relativistic QCD framework [18]. 
The COM allows the formation of a charmonium from a color-octet cc pair 
with one or some soft gluon emissions. Figure 4 shows an example of the 
lowest order production of a, J/ip with the COM from the gluon fusion sub- 
process. The cc pair, which is in ^5*0 (or ^Pj) and the color-octet state, forms 
a J/i/j with a soft gluon emission. Using appropriate color-octet matrix el- 
ements, which are additional free parameters needed to be extracted from 
experimental data, the COM has successfully reproduced pr distributions 
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Figure 3: An example of the lowest order diagram for direct J/ip production 
from gluon fusion with the color-singlet model. Incoming protons are omitted 
in the figure. The cc pair is in the color-singlet state. 

at CDF [19-21] as shown in Fig. 5 and total cross sections at lower-energy 
experiments [22-24]. COM predictions for the relative yields for each char- 
monium state (for example, Xc to direct J/ip ratio) arc also consistent with 
the experimental data [23], which will be discussed in Appendix A. However 
extraction of these matrix elements is still controversial and therefore large 
ambiguities are left for the prediction with the COM. 




Figure 4: An example of the lowest order diagram for direct J/ijj production 
from gluon fusion with the color-octet model. Incoming protons are omitted 
in the figure. The cc pair is in the color-octet state. 

Hadro-production data on J/ ip and other charmonia are available in both 
fixed-target experiments and collider experiments in a wide energy range 
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Figure 5: Transverse momentum differential cross sections for prompt (ex- 
cluding 6-quark decays) J/ip production measured by the CDF experiment 
compared with theoretical predictions. The dashed curve depicts the direct 
color-singlet contribution. The dot-dashed curve illustrates the contribution 
of the ^Si octet state of the cc pairs and the dotted curve denotes the com- 
bined contribution of the ^5*0 and ^Pj octet states. The sohd curve equals 
the sum of the color-singlet and color-octet contributions. All curves are 
multiplied by the muon branching ratio B{J/ip — > 
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(6.1 GeV < < 1.8 TeV) [25-48]. Measurements of such as total and pr- 
differential cross-sections, polarization, and relative yields of each charmo- 
nium have promoted better understanding of production mechanism. How- 
ever more data on different observables at different energies are required 
since currently none of the theoretical models can successfully explain all the 
experimental data at all energies. 



Our measurements of J/ip in p+p collisions using the PHENIX detector 
at RHIC will result in another critical input for understanding of hadro- 
production of J/ijj. The energy region of RHIC {^/s — 200 to 500 GeV) 
has never been exploited yet by previous experiments. Thanks to excel- 
lent lepton-identification capabilities of two independent spectrometers of 
PHENIX which cover both the central rapidity region {\y\ <0.35) and the 
forward region (1.2 < y <1.2), we are able to extract total cross-section for 
J/ ip production with less ambiguity in the extrapolation to unmeasured kine- 
matical region, unlike higher-energy {y/s > 630 GeV) collider experiments 
with a pt cut (pr > 4 or 5 GeV/c) and limited kinematical coverage in the 
forward region. Therefore, we can obtain total cross-section results at the 
highest energies, which are sensitive to the production models. 

In addition, our measurements of differential cross sections will provide 
further information to understand perturbative QCD. Our wider kinemati- 
cal coverage enables us to measure rapidity-differential cross section which is 
reflected by gluon distribution function g{x) in the proton. Transverse mo- 
mentum (pt) distribution at lower px, which dominates our inclusive yield, 
and the average value of pt are sensitive to intrinsic transverse-momentum 
of partons (kx) ■ 

RHIC has started a physics run from year 2000. In the 2001-2002 run 
period (Run-2), both Au+Au collision and p-|-p collision data have been ac- 
cumulated, where J/ particles have been successfully detected via both the 
e+e" and decay channels. In this paper, measurements of total and 

differential cross sections for inclusive J/ip production in p+p collisions at 
^/s — 200 GeV are presented. In section 2 and 3, details of the experimen- 
tal setup and analysis procedure are described respectively, focusing on the 
fi~^fi~ channel measurement. In section 4, results arc discussed together with 
theoretical predictions as well as results of other experiments, followed by 
the conclusion. 



R D. Sato 



13 



2 Experimental setup and data acquisition 

2.1 The RHIC accelerator complex 

The Relativistic Heavy-Ion Collider (RHIC), located in Brookhaven Na- 
tional Laboratory, Upton, New York, is capable of accelerating a wide variety 
of nuclei and ions from protons to Au (gold) nuclei up to 250- GeV energy 
for protons (or 100-GeV per nucleon for Au) using two independent rings 
and colliding them at six interaction points. The design luminosities are 2 x 
10^^ cm^^ scc^^ for Au beams and 2 x 10^^ cm~^ sec~^ for proton beams (2 
X 10^^ cm~^ sec~^ in an enhanced mode) at the top energy. 

RHIC also has a capability of accelerating polarized protons. In the 
Run-2 period (year 2001-2002), RHIC has successfully achieved transversely- 
polarized proton-proton coUisions at = 200 GeV and delivered an inte- 
grated luminosity of roughly 1 pb~^ to each experiment. In this paper we 
will present unpolarized cross sections for the J/i(j production, that is, spin- 
averaged results. Therefore, descriptions of the polarized proton acceleration 
are not described here. 



12 o'clock 




Figure 6: The proton acceleration to RHIC 

Attached figure 1 shows an aerial view of the RHIC accelerator complex. 
The route for the proton acceleration to RHIC is shown in Fig. 6. Starting 
from the polarized ion source, (polarized) protons are accelerated through 
Linac, Booster and AGS then injected into both rings of RHIC, whose cir- 
cumference is 3.834 km. The Blue ring runs clockwise and the Yellow ring 
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runs counter-clockwise. There are currently 120 bunch buckets in each ring 
whose interval is 106 nsec (or 9.4 MHz frequency). In Run-2, only about a 
half of them (54 bunch buckets) have been filled. Typical bunch length was 
2 nsec (60 cm) in rms during the p+p run. 

Experiments are located at the interaction points in RHIC where two 
bunches in each ring coUide at an angle of 0°. There are six interaction points 
called 12, 2, 4, 6, 8, and 10 o'clock respectively starting from the north and go- 
ing clockwise. Table 1 shows the experiments at each interaction point. Each 
experiment demonstrates its unique feature. STAR and PHENIX are the 
largest experiments at RHIC each with more than 400 collaborators. STAR 
(Solenoid Tracker At RHIC) tracks and identifies charged particles with a 
time projection chamber covering a large solid angle. PHENIX (Pioneering 
High Energy Nuclear and Ion cXperiment) was designed to measure hadrons, 
leptons and photons in both high multiplicity and high rate environments. 
BRAHMS measures hadrons over wide ranges of rapidity and momentum 
using two magnetic spectrometers. PHOBOS consists of a large number of 
silicon detectors surrounding the interaction region to measure charged par- 
ticle multiplicities even in the most central Au+Au collisions. The pp2pp 
experiment is aiming at measuring p+p total and elastic cross sections. 



Interaction Point 


Experiment 


(o'clock) 




12 




2 


pp2pp, BRAHMS 


4 




6 


STAR 


8 


PHENIX 


10 


PHOBOS 



Table 1: Experiments at RHIC 



2.2 Coordinates and formulae 

In this subsection, coordinate system and formulae used to describe the 
experimental setup and results are introduced. 

Coordinate system at PHENIX 

The z axis stands for the beam line which runs straight in the experi- 
mental area, where positive z points to the north. The polar angle 6 and 
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azimuthal angle are defined with respect to the z axis. The north (south) 
direction is usually defined as ^ = (180) degrees and the west (east) direc- 
tion is defined as (p — (180) degrees. 

Rapidity 

Rapidity y of a particle is defined as 

1, E + p, 

y = - In — 

2 E-p, 

where E is the energy and Pz is the ^-component of momentum of the particle. 

Pseudo rapidity 

Pseudo rapidity 77 of a particle is defined as 

_ 1 Ptot +Pz . . 

r] = - In = — in tan - 

2 Ptot -Pz 2 

where ptot is the scalar value of momentum and 9 is the polar angle of the 
particle direction. In the massless limit, y reaches rj. 

In the definition above, y(r)) of a particle going to the south is negative. 
However, absolute (positive) values of y{r]) will be often used for those par- 
ticles in this paper, since results with the spectrometer on the south side 
of PHENIX will be mainly described and physics should be symmetric with 
respect to y{ri) = 0. 

2.3 The PHENIX experiment overview 

PHENIX is one of the largest experiments at RHIC, located at the 8- 
o'clock interaction region. Attached figure 2 shows a schematic view of the 
PHENIX experiment. PHENIX was designed to measure leptons, photons 
and hadrons in both high-multiplicity heavy-ion collisions and high event- 
rate p-|-p collisions. 

There are two independent spectrometers in PHENIX which cover differ- 
ent pseudo-rapidity regions. Two Central Arms, East and West Arms, cover 
the pseudo-rapidity range of \ri\ < 0.35 with a quarter azimuth for each Arm 
and measure electrons, photons and hadrons. Two Muon Arms, North and 
South Arms, cover 1.2 < 77 < 2.4 and —2.2 < 77 < —1.2 respectively with a 
full azimuth and measures muons. 

There are three magnets in PHENIX. The Central Magnet provides an 
axial magnetic field for the Central Arms while two Muon Magnets produce 
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a radial field for each Muon Arm. Attached figure 3 shows magnetic field 
lines inside the magnets. A part of the Central Magnet steel and a copper 
spacer mounted on each side of it (called a copper nosecone) work as a hadron 
absorber of about five interaction-length for the Muon Arms. The positions 
of the nosecones, z = ± AO cm, determine a useful vertex region for physics 
events (|2;| < 30 to 40 cm depending on analysis). 

The Central Arms consist of three kinds of tracking chambers (Drift 
Chambers, Pad Chambers and Time- Expansion Chambers), Ring-Imaging 
Cerenkov Counters for the electron identification, Time-Of-Flight detectors 
for the particle-identification for hadrons, and Electro-Magnetic Calorime- 
ters for measuring energies of electrons and photons. These detectors are 
positioned radially with respect to the 2;-axis extending from 2 m to 5 m. 
Details of operation and performance of the Central Arm detectors can be 
found in [49-51]. Descriptions of the analysis of the J/ijj ^ e^e~ channel 
using the Central Arms are found in [52]. 

In addition to these four Arms, there are three kinds of counters to trigger 
p-|-p interactions which will be described in the next subsection followed by 
the subsection dedicated to description of the Muon Arms. 

2.4 Interaction trigger counters 

Three kinds of interaction-trigger counters (ITC) have been used to trig- 
ger p-|-p inelastic events and find vertices during the Run-2 pH-p period. They 
are Beam-Beam Counters (BBC), Normalization Trigger Counters (NTC) 
and Zero-Degree Calorimeters (ZDC). Table 2 summarizes acceptance and 
performance of each ITC. Since trigger efficiency of the ZDC for p-|-p inelas- 
tic events is small (0.01), it was not used for the J/ip analysis, hence will not 
be described here. Detailed information on the ZDC is found in [53] 



Counters 


Acceptance 


Sensitive particle 
type 


Typical 
efficiency 


Typical vertex 
resolution (cm) 


BBC 


3.0< \r]\ <3.9 
13 > 0.7 


charged particles 


0.5 


2 


NTC 


1.1< \r]\ <2.8 


charged particles 


0.6 


10 


ZDC 


rj >6.2 


neutral particles 


0.01 


10 



Table 2: Acceptances and performances of the interaction-trigger counters. 
Efficiencies and resolutions are for p-|-p inelastic events. 
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2.4.1 The Beam-Beam Counters 

Two Beam-Beam Counters (BBC) have been used for the primary trigger 
and vertex counter for the p+p interactions. They are positioned 1.4 me- 
ters away from the interaction point along the beam axis on each side and 
cover from 2.4 to 5.7 degrees (3.0 < \r]\ < 3.9) with a full azimuth. They 
determine the vertex position of an event from the time difference between 
two counters. Attached figure 4 shows a photograph of a BBC. Each BBC 
consists of 64 hexagonal Cerenkov radiators each of which is mounted on a 1" 
photomultipher tube (PMT). They are sensitive to charged particles whose 
f3 = v/c is greater than 0.7. 

With a beam test, intrinsic timing resolution of 50 psec was obtained for 
one module [54]. In Au+Au collisions, better vertex resolution is expected 
because hit multiplicity is higher. Actually 0.5-cm vertex resolution has been 
obtained in Run-2 Au-|-Au collisions. In p-|-p coUisions, about 2-cm vertex 
resolution is expected which is still good enough not to worsen invariant mass 
resolution for J/ip particles measured in the Muon Arm. Trigger efiiciency 
for p+p inelastic events is expected to be about 0.5, which will be discussed 
in section 3. 

Beam-Beam Local-Level- 1 

Hit time information of all PMTs is sent to the Beam-Beam Local-Level- 1 
(BBLLl) board where a trigger decision is made. Online vertex position is 
obtained from hit time information without pulse-height corrections (threw- 
ing corrections). During the p-|-p run, a trigger is fired when there is at 
least one hit on both sides of the BBC counters and online 2;- vertex position 
\zvtx\ < 75 cm which is sufficiently large compared to offline vertex cuts (30 
to 40 cm). Online vertices (without slewing corrections) and offline vertices 
(with slewing corrections) agree within the accuracy of the vertex determi- 
nation (2 cm). 

2.4.2 The Normalization Trigger Counters 

Normalization Trigger Counters (NTC) have been introduced to increase 
trigger efficiency for p-|-p inelastic (including diffractive) events which is 
about 0.5 with the BBC only. 

Each NTC is located on top of each side of the noseconc which is 40-cm 
away from the interaction point. Pseudo-rapidity coverage is 1.1 < \ri\ < 2.8. 
It consists of four fan-shaped scintillators (called quadrants) each of which 
is mounted on a PMT to collect scintillation light emitted when a charged 
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particle traverses. Single-particle detection efficiency obtained with a beam 
test is about 90% [55]. 

A simple NIM-logic makes NTC triggers using hit time information of all 
quadrants, which are sent to the Global-Level- 1 board (see section 2.7). A 
trigger is fired when at least one quadrant on both sides has a hit. Trigger 
efficiency for p-|-p inelastic events is estimated with a simulation to be 60% 
by the NTC itself and go up to 74% when combined with the BBC where 
statistical errors of the simulation is 1%. Although NTC-triggered events 
without a BBC-trigger have not been used to increase the number of J/ip^s 
because of its poor vertex-resolution (10 cm), they are used to confirm trigger 
efficiency of the BBC which will be discussed in section 3. 

2.5 The Muon Arms 

The PHENIX Muon Arms were designed to detect muon pairs from decays 
of vector mesons and bosons produced in the forward rapidity region 
(-2.2 < y < -1.2 for the South Muon Arm and 1.2 < y < 2.4 for the North 
Muon Arm) as well as single muons from (semi-)leptonic decays of open 
heavy flavors and W bosons. Attached figure 5 shows a schematic view of a 
Muon Arm. Since only the South Arm was operational in Run-2, following 
descriptions specifically focus on it. 

In this forward region, hadron background is relatively larger than that 
in the central region, since hadrons can be produced in soft processes with 
wider rapidity distribution than hard processes, while signal muons such as 
from heavy flavors are produced in hard processes. Good hadron rejection 
while keeping good momentum resolution is achieved with the following com- 
ponents of the Muon Arms and PHENIX. 

1. Pre-rejection of hadrons with a five- interaction- length absorber of the 
nosecone and Central Magnet. Interaction length was determined not 
to degrade momentum resolution for low momentum (2 to 5 GeV/c) 
muons and to keep good mass resolution for J/'^'s. Hadron rejection 
factor of about 100 is achieved here. 

2. Measurement of particle momentum with a magnetic spectrometer (Muon 
Tracking Chamber inside the Muon Magnet) 

3. Further rejection of hadrons with an array of coarse-segmented tracking 
chambers and absorbers (Muon Identifier). Another factor of about 30 
is achieved for the hadron rejection. 
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Figure 7 shows integrated nuclear interaction-length (Xint) in the South 
Muon Arm as a function of the distance from the interaction point in the z 
direction. At the last MuID gap (gap 5), Xint becomes 9.65. The minimum 
Pz (z-component of momentum) for a muon to reach gap 5 is 2.5 GeV/c. 
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Figure 7: Integrated nuclear interaction-length of the absorbers in the South 
Muon Arm as a function of the distance from the interaction point in the z 
direction. Vertical lines indicate rough positions of the chambers. Hatched 
areas represent absorber materials. 



2.5.1 The Muon Tracker 

The PHENIX Muon Tracker (MuTr) comprises three stations of tracking 
chambers inside the Muon Magnet as shown in Fig. ??. Its design was driven 
by requirements from both heavy-ion physics and spin physics. The separa- 
tion of each charmonium or bottomonium state from the others, J/ip from ip^ 
for example, is essential to find a QGP signal, since the degrees of suppression 
for each state are expected to vary because of different binding radii. For 
spin physics, charges of high-p^ {pt > 20 GeV/c) muons from W and Z bo- 
son decays are needed to be identified [8]. To satisfy the requirements above, 
100-/im resolution is needed for chamber resolution. In addition, multiple 
cathode-strip orientations and read-out planes are required for each station 
to reconstruct tracks efficiently even in the most central Au-|-Au events. A 
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MuTr electronics design was also driven by the requirement of 100-/im reso- 
lution measurements. 

With 100-/im position resolution, momentum resolution Ap/p = 3 to 5% 
is achieved for 2 to 10 GeV/c muons as shown in Fig. 8. For low momentum 
(2 to 5 GeV/c) muons, multiple scattering is the dominant factor to smear 
muon momenta, whereas position resolution of chambers becomes dominant 
for high momentum (above 10 GeV/c) muons. Polar angle dependence of 
momentum resolution for high momentum muons is due to the difference in 
magnitudes of the magnetic field inside the South Muon Magnet as shown in 
Table 3. 




Figure 8: Momentum resolution for muons with different polar angles in 
the South Arm acceptance as a function of pz obtained with a simulation. 
Typical statistical errors of the simulation is 5%. 



Polar an^ 


s;le in degrees 


165 


160 


155 


150 


Integrated ma} 


2;netic field in T-m 


0.774 


0.494 


0.344 


0.255 



Table 3: Integrated magnetic field inside the South Muon Magnets. 



Mechanical design 

All chamber planes inside the Muon Magnet are perpendicular to the z- 
axis. Their ^-positions for each station are shown in Fig. 7. The magnetic 
field is in the radial direction so that charged particles from the interaction 
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vertex bend primarily in the azimuthal direction, which is perpendicular to 
the direction of cathode strips to determine the track positions. Typical 
bend from the straight line at station 2 is 1 cm for medium-momentum (5 
to 10 GcV/c) muons. 

In one MuTr station there are three (or two) gaps each of which con- 
sists of two cathode-strip planes and one anode- wire plane in between with a 
3.2-mm anode-cathode spacing. Attached figure 6 shows a cross section of a 
MuTr station. Station 3 has only two gaps because it can exploit additional 
two-dimensional position information of MulD roads. All gaps are divided 
into octants electrically as shown in attached figure 7, together with octant 
numbers. For station 2 and 3, octants are also the unit of mechanical assem- 
blies while quadrants are for station 1. Attached figure 8 shows a photograph 
of a station-2 octant. An octant is further divided into two half-octants in 
the middle, in each of which directions of anode wires and cathode strips are 
fixed. 

A cathode plane consists of 5-mm width strips with alternate readout to 
avoid cross-talks between them. An anode plane is an alternating structure 
of 20-;um gold-plated tungsten sense wires and 75-fim. gold-plated Cu-Be field 
wires with a sense wire spacing of 10 mm. Anode wires run in the azimuthal 
directions while cathode strips run in the radial directions. The direction 
of cathode strips in one plane of each gap is perpendicular to that of the 
anode- wires (called a non-stereo angle plane), with a goal of 100 /xm position 
resolution. The direction of cathode strips in the other plane of the gap 
(a stereo angle plane) is tilted by 3.25 to 11.25 degrees depending on gap 
and station, which are summarized in Table 4 and illustrated in attached 
figure 9 for station- 1 planes. Resolutions of stereo-angle planes are worse 
(300 /im) but they are needed to determine 3-D positions of hits and reject 
ghost tracks. 

Specific technologies were used for each station to produce a cathode 
pattern to an accuracy of better than 25 /im; photolithography for station 1, 
electro-mechanical etching for station 2 and mechanical routing for station 3 
[56] . A unique wire laying apparatus was designed and implemented for each 
station. 

In order to reduce the multiple scattering in the spectrometer which de- 
grades momentum resolution, thickness at the station 2 detector was required 
to be less than 10~^ of a radiation length. To meet this requirement, the 
station-2 octant cathodes were made of etched 25-//m copper coated mylar 
foils. The thickness of the copper coat is 600 A. As a result, the total thickness 
of 8.5 X 10"^ radiation lengths is achieved which satisfies the requirement. 

The chamber gas mixture was 50% Ar + 30% CO2 + 20% CF4 with a gas 
recirculation system included in normal operation. The nominal high-voltage 
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Station 


Gaj) 


angle (degree) 




1 


-11.25 


1 


2 


+6 




3 


+ 11.25 




1 


+7.5 


2 


2 


+3.75 




3 


+11.25 


3 


1 


-11.25 




2 


-11.25 



Table 4: Relative angles of the cathode strips in stereo-angle planes with 
respect to the non-stereo angle (radial direction). Positive signs represent 
the positive 4> direction (counter-clockwise) for all half-octants in station 1 
and half-octant in station 2 and 3, while the negative (f) direction for the 
others. 



potential applied to anode wires was 1850 V with a gain of approximately 
2 X 10^ 

To maintain good momentum resolution, an optical alignment system has 
been installed to calibrate initial placement of the chambers, and to monitor 
displacement of the chambers during their operation to ±25 /xm. There 
are seven optical beams surrounding each octant chamber, consisting of an 
optical- fiber light source at station 1, a convex lens at station 2 and a CCD 
camera at station 3 which are shown in attached figure 10. 

Electronics design 

Attached figure 11 shows a schematic diagram for the MuTr Front End 
Electronics (FEE). Raw chamber signals are continuously amplified by CPAs 
(Charge Pre- Amps) and stored in AMUs (Analog Memory Units) with 64- 
event buffers with the 10-MHz beam clock. Upon receipt of a level-1 trigger 
bit from a Granule Timing Module (GTM), stored samples of all channels 
are digitized by 11-bit ADCs (Analog to Digital Converters) and the results 
are sent to a Data Collection Module (DCM) through a Front End Module 
(FEM). GTMs and DCMs are described in section 2.6. 

Four ADC samples are used to determine the amount of the charge de- 
posited on a strip to reject noise hits as much as possible. Each sample is 
measured for the duration of 100 nsec. The second sampling starts 400 nsec 
after the first sampling ends. Second to fourth samplings are consecutive. 
Timing has been set so that the third sample comes to the peak of the pulse 
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as shown in Fig. 9. Relative charges (or ADC counts) of these samples have 
been monitored online to guarantee peak positions not to move around from 
the third sampling. The amount of the peak charge of a strip is obtained 
offline as the average of the second to fourth samples with a pedestal sub- 
tracted. These four samples are converted within 40 usee per event. 
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Figure 9: Typical ADC counts for each sampling together with a typical 
signal pulse (dotted line). The horizontal axis is for the relative sampling 
time. One ADC count corresponds to about 0.5 fC (10~^^ Coulomb). 

Strip by strip calibration is crucial for good position resolution, since 
a position of a muon track is determined by fitting charges on typically 
2 or 3 consecutive strips induced by the track, which will be described in 
section 3.2.2. A calibration system has been implemented to inject pulses 
into all of the chambers. Four wires in each chamber gap, which span the 
entire width of the cathode planes, are sent a square pulse from a digital to 
analog converter (DAC), thus inducing a charge on all cathode strips in a 
given gap simultaneously. Several different pulse amplitudes are sent to the 
chambers and many events are collected at each amplitude so that relative 
gains of the cathode strips can be determined over the entire range of the 
electronics. Pedestals are monitored by collecting calibration data with the 
DAC amplitude set to zero. 

To meet the design requirement of 100-/im resolution, rms noise at the 
input to the preamps is required to be 0.5 fC for a typical pulse of 80 fC, 
which is achieved in a test bench measurement described next. 
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Integrated performance 

Integrated performance of chambers and electronics has been studied in 
a cosmic-ray test in a test stand with one station-2 chamber and its full 
complement of electronics, and in readout of the entire South MuTr system 
prior to Run-2. The cosmic-ray test data showed that the system was ca- 
pable of meeting the noise specifications and that 100-/im resolution could 
be achieved. The noise specifications have been met on the full South MuTr 
system and the system has been shown to be robust over several months of 
data taking. 

The cosmic-ray test was performed with one station-2 chamber, 960 chan- 
nels of production front-end electronics, the same high-voltage and low- 
voltage distribution system that is used in the final system, and with a copy 
of the PHENIX data acquisition system. The noise specifications of 0.5 fC 
(1 ADC count) were met, as can be seen in Fig. 10, where the rms values of 
the pedestals on all readout channels are shown. Noise environment in situ 
has turned out to be similar to this test except for chambers in station 2, 
which has caused degradation of the mass resolution for J/ip in Run-2. 




Figure 10: Measurement of rms noises for 128 typical channels in units of 
ADC counts (1 ADC count is roughly 0.5 fC) obtained in the cosmic-ray 
tests for a station-2 octant. 



Two scintillators, one on either side of the station-2 chamber, were used 
to provide a trigger for cosmic rays going through the chamber. The data 
collected from this trigger were searched for clusters in each cathode readout 
plane, the clusters were fit to extract the centroid strip positions, and 5 out 
of 6 readout planes were fit to a straight line and projected to the sixth. 
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central non-stereo readout plane. A cut was placed on the straight line fit 
to only the select tracks which were approximately perpendicular to the face 
of the chamber and the difference between the projected straight-line fit and 
the measured position on the sixth plane was plotted. The result is shown in 
Attached figure 12, where a resolution of approximately 100 fim was achieved 
when the projection error, position resolution for a track obtained with the 
other five planes, is removed from the residual. 

Each individual channel's gain, pedestal and variation (or noise) in the 
pedestal were measured. The dynamic range in the charge measurement 
of the system was verified and long runs demonstrated the stability of the 
optical links from FEMs to a DCM. Figure 11 shows the residuals from a gain 
measurement during the commissioning period. The residuals are shown to 
be consistent with a linear gain to a few ADC counts over the operable ADC 
range. 
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Figure 11: Residuals (in ADC counts) from a straight-line fit to the charge 
measured on a given strip of the ADC vs the DAC pulse amplitude applied 
to the calibration wires. Two different strips from a given cathode plane 
are shown. Approximately one half of the ADC range (of 2047 channels) is 
shown and is linear to a few ADC counts and constant over time. 



2.5.2 The Muon Identifier 

The Muon IDentifier (MuID) consists of five layers of chambers inter- 
leaved with steel absorbers. Each chamber plane is called gap 1 through 



26 



J/V' Production in p+p Collisions at ^/s = 200 GeV 



gap 5 starting from the nearest gap to the interaction point. The MuID is 
used for separating muons from charged hadrons and other background as 
well as providing triggers for single muons and dimuons (muon pairs). It 
gives another hadron rejection factor of about 30 in addition to the Central 
Magnet and nosecone (about 100), thus reducing the mis-identification rate 
for the punch-through hadrons to 3 x 10~^. This is much smaller compared 
to the irreducible hadron weak-decay background (tt — > /ii' and K fiv) in 
fiight before the nosecone for near-threshold momentum (3 GeV/c) muons 
(3 X 10~^). Minimum for a muon produced at the interaction vertex to 
reach the MuID is 1.8 GeV/c and 2.5 GeV/c to penetrate through it. 

Segmentation of the absorber into multiple layers improves the measure- 
ment of the trajectory in the MuID gaps (chamber layers) for low momentum 
muons, which is desirable to increase the acceptance for the (f) meson detec- 
tion. The segmentation chosen is 20 cm, 10 cm, 10 cm, 20 cm and 20 cm 
starting from the South Magnet backplate as shown in Fig. 7. In each gap 
between these absorbers, chamber panels are installed. 

Mechanical Design 

One MuID plane, or a gap, consists of six panel structures (called MuID 
panels) as shown in Fig. 12, into which chambers are assembled in both 
horizontal and vertical orientations. The upper figure of Fig. 13 shows the 
cross section of a MuID panel. 

larocci-type plastic tubes have been chosen as MuID chambers for longevity, 
robustness and low cost to cover a large area (13 x 10 m^ for each Arm). One 
tube has eight cells with a 9 x 9-mm^ cross section each of which has a gold- 
coated CuBe anode-wire with a 100-/im diameter at its center. The lower 
figure of Fig. 13 shows the cross section of an larocci tube (in this figure, two 
tubes, or a two-pack is shown which will be explained later). The cathode 
wall is made of polyvinyl chloride (PVC) coated with carbon. Length of a 
tube varies from 2.5 m to 5.6 m depending on its position to be installed. 

Tubes are operated in the proportional mode with 4300 to 4500 V poten- 
tials and isobutane-C02 mixed gas, where a gain of approximately 2 x 10^ is 
achieved. As test-bench results, 92 ± 1 % efficiency, which was also measured 
in the beam test to be described later, and 80 usee (nano-seconds) drift time 
width were obtained. To achieve a better efficiency and faster drift time, 
two tubes with a half cell (5-mm) shift consist one channel (called a 2-pack) 
as shown in Fig. 13. For a 2-pack, 97 ± 1 % efficiency and 60 nsec (nano- 
seconds) drift time have been obtained [57]. Drift time is faster enough than 
the bunch-crossing interval (106 nsec), so that a level-1 trigger can uniquely 
determine the bunch crossing in which an event has occurred. Gas mixture 
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Figure 12: Configuration of MuID panels in one gap. The beam axis runs in 
the middle perpendicularly to the paper. 

ratio is adjustable between to 25% for isobutane. In Run-2 it has been set 
to 7% to meet the non-flammable requirement of PHENIX. 

larocci tubes are assembled in a panel structure called a MuID panel 
which consists of Al frames, cover plates and a mid-plate (shown in the upper 
figure of Fig. 13). larocci tubes are glued to both sides of a 3-mm width Al 
mid-plate with double-sided tapes at 8.4-cm intervals for both horizontal and 
vertical orientations. Each orientation has two chamber layers shifted by a 
half cell (5 mm) to make 2-packs. Those layers have independent chains for 
gas and high-voltage supplies to minimize the number of dead channels. 

There are two kinds of the panel size which are 5.6 x 5.2 m^ (large panel) 
and 4.4 (or 4.2) x 2.9 m^ (small panel). In one gap, four large panels and 
two small panels are installed as shown in Fig. 12. Adjacent panels overlap 
with each other so that there is no inactive area between them. 

Electronics Design 

A passive OR of signals of 16 wires in a 2-pack is read out and amphfied 
by a factor of 150 with an in-panel amplifier on a high-voltage distribution 
board which also provides high voltage to each tube. Signals are then sent 
to the MuID Front-End Electronics (FEE) where they are again amplified 
by a factor of 3, discriminated and stored in the data buffer. Discriminator 
threshold values have been set to 90 mV to minimize the number of noise 
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plastic tube 




side frame mid-plane 

Cross section of the MulD panel 




83.5mm 



Cross section of the plastic tube(2-pack) 

Figure 13: Cross section of a MuID panel (upper) and an larocci tube (lower). 
In one panel-orientation, there are two layers with independent gas and high- 
voltage chains. Two adjacent tubes in those layers make one signal channel 
to be read out. 
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hits while keeping good efficiency (typical pulse height of signals is 500 mV 
to 1 V). Upon receipt of a level-1 trigger, all digitized bits are sent to a DCM. 
There is additional output of the MulD FEE called "pseudo-trigger output" 
which is grand logical OR of a certain fraction of channels. They are used 
for the NIM-logic level-1 trigger which is described in section 2.5.3. 

Integrated performance 

To confirm hadron-rejcction performance and muon detection efficiency 
of the entire MulD system experimentally, a beam test was performed at 
KEK-PS ^ using the same type of larocci tubes and steel absorbers in the 
same configuration as in PHENIX. 

Attached figure 13 shows the experimental setup for the beam test. Pion 
and muon beams were produced by bombering an inner target with 12-GeV 
proton beams, then 1-4 GeV/c momentum was selected with the magnet 
system and delivered to the experimental area. Four scintillation counters 
(STl to ST4) defined the beam. Three gaseous Cerenkov counters (GCl 
to GC3) identified pions and muons. GCl and GC2 were pressured to dis- 
tinguish between pions and muons and GC3 to distinguish between muons 
and electrons. Beam qualities obtained were better than 99% for muons and 
better than 99.9 % for pions excluding their weak decays (tt^ — > /x^t^) after 
GC2. Five iron slabs were used with a width of 10 cm, 10 cm, 10 cm, 20 cm 
and 20-cm respectively staring from the upper stream. An additional 10-cm 
(20-cm) plate was added in front of the first layer to simulate the backplate 
of the South (North) Muon Magnet. The numbers of larocci tubes used were 
3, 3, 5, 7 and 9 starting from the first gap for each orientation. They were 
optimized to 3a dispersion of 2-GeV/c momentum muons due to multiple 
scattering. 

As a result, muon detection efficiency of 86 ± 2 % has been obtained for 
1.8-2.5 GeV/c muons with a small momentum dependence. This is slightly 
lower than the test bench result described before (92 ± 1 %), which is ex- 
plained by additional inactive volumes between tubes. Attached figure 14 
shows the results of pion mis-identification rate as a function of pion mo- 
mentum, which is consistent with a GEANT [58] simulation including weak 
decays into muons. For the South Arm, mis-identification rate for 4-GeV/c 
momentum pions (about 5 GeV/c at the interaction region) has been deter- 
mined to be 0.04 excluding decays. Multiplying the rejection factor of the 
nosecone and Gentral Magnet {e~^ ~ 7 x 10~^), the net mis- identification 
rate of 3 x 10~^ is obtained, which satisfies the design value. 



^Proton Synchrotron at High Energy Accelerator Research Organization, Japan 
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2.5.3 The MuID NIM-logic trigger for LVL-1 

A coarse-segmented trigger system called the MuID NIM-logic trigger^ 
was used to trigger muons in the South Muon Arm during the Run- 2 p-|-p 
period. Attached figure 15 shows a schematic view of the system. 

As shown in the figure, all the MuID planes are divided into four regions 
(called quadrants) by horizontal and vertical lines through the middle. A 
trigger decision for a quadrant, that is, whether a muon has passed or not, 
is made by taking the coincidence of fired planes of the quadrant. Since the 
number of input channels of the trigger circuit is limited, only four gaps out 
of five were used. Deeper gaps (gap 3,4 and 5) were included since they are 
more important for muon identification than shallower gaps. Gap 1 was also 
included since (1) better tracking performance due to longer tracking length 
and (2) better rejection of cosmic rays and very low angle (^ < 9 degrees) 
particles. The chance for a muon to cross more than one quadrant is small. 
Each quadrant is further divided into two sectors for both horizontal and 
vertical tube orientations. Figure 14 shows the segmentation for each gap. 
An output of each sector, grand OR of all channels inside, corresponds to a 
pseudo-trigger output of a ROC. in total, 16 pseudo-trigger output signals 
are sent to the trigger algorithm. 
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Figure 14: Segmentation of a MuID plane for the NIM-logic trigger. Bold 
lines divide quadrants and thin lines divide segments in each quadrant. Cap- 
ital letters represent positions of each segment and used in the text. 

Trigger decisions are made by LeCroy 2372 Memory Lookup Units (MLUs) 
for each quadrant. Two MLUs arc prepared for each quadrant for both 
"deep" and "shallow" triggers. A deep trigger requires hits in all four gaps 
used and a shallow trigger requires hits up to gap 3. Hit patterns are required 
to point to the event vertex. Hit patterns such as A-A-B-B and A-C-D-D are 
accepted but B-B-A-A and D-A-A-C are not, for example, where four letters 
represent hit segments in a quadrant shown in Fig. 14 for each gap starting 

^So named because the trigger circuit is constructed using NIM- and CAMAC-standard 
electronics. 
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from the first gap. An example of an accepted-pattern is shown in attached 
figure 16. To minimize the loss of trigger efficiency due to finite chamber 
efficiencies, the algorithm allows some gaps to miss hits. For a deep trigger, 
6 out of 8 gaps (including both orientations) are required to have a hit and 3 
out of 4 gaps for a shallow trigger. Extra hits are allowed, thus no efficiency 
loss is expected due to background hits. 

Total 8 quadrant trigger signals are sent to another MLU, called the 
"Decision MLU" , which counts the number of triggered quadrants and issues 
five trigger signals: 



(1) single shallow (IS) - one shallow quadrant is fired, 

(2) double shallow (2S) - two or more shallow quadrants, 

(3) single deep (ID) - one deep quadrant, 

(4) deep-shallow (ID IS) - one deep quadrant and one or more shal- 

low quadrants, and 

(5) double deep (2D) - two or more deep quadrants. 



which are sent to the Global Level- 1. A shallow quadrant-trigger is fired 
always when the deep quadrant-trigger is fired for that quadrant. Therefore, 
the following relations hold true: 2D C IDIS C ID C IS where A C B 
denotes that the trigger B is always fired when the trigger A is fired. The 
reasons to have both shallow and deep triggers for each quadrant are (i) to 
achieve higher J/ip ^ fi~^fi~ efficiency obtained with the IDIS trigger than 
the 2D trigger by 30% and (ii) to study hadron punch-through background 
with the IS and 2S triggers (not described in this paper). 



2.6 The PHENIX DAQ system 

Figure 15 shows a schematic view of the PHENIX Data Acquisition 
(DAQ) System. All the subsystems are equipped with the timing mod- 
ules called Granule Timing Modules (GTMs). The Master Timing JVIodule 
(MTM) dehvers the 9.4-MHz RHIG clock to GTMs. When a trigger is issued 
by the Global Level-1, the clock and trigger bit are transported over optical 
fibers to the electronics of each detector, the Front End Module (FEM). The 
detector records the data in raw digitized format and transports the data 
packets over optical fibers to the Data GoUection Module (DCM). The data 
are recorded in a buffer disk in the PHENIX control room where data quality 
is monitored online. At the maximum 60 mega-bytes per second has been 
achieved as a rate of recording data in the disk which corresponds to about 
1200 events per second for p-hp data whose size is about 50 kilo-bytes per 
event. The data are then transferred and stored in a huge storage tape with 
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1.2x10^^ bytes capacity (sharing with other RHIC experiments) for offline 
computing such as cahbration and track reconstruction. In Run-2 PHENIX 
has recorded about 2x10^'^ bytes for the p+p data. 




Figure 15: The PHENIX data acquisition system, triggered by an event in 
the Zero-Degree Calorimeters (Left in diagram) 



2.7 Global Level-1 triggers 

The level-1 triggers of PHENIX are issued by the Global Level-1 (GLl) 
module, for which local level-1 and NIM-logic triggers from various subsys- 
tems are input. The maximum DAQ rate (about 1 kHz) and the number of 
triggers (10 or more) limit the rates for each GLl trigger to about 100 Hz or 
less. Level-2 triggers were not used during the p-|-p Run. 

2.7.1 Minimum bias trigger 

Logical OR of the BBC (BBLLl) and NTC triggers, called the "Minimum 
Bias (MB)" trigger was used to trigger p+p inelastic events with about 70% 
efficiency. MB-triggered events are useful for studying detector efficiency 
and background as well as extracting physics. At the highest luminosity of 
RHIC achieved in Run-2 (10^° cm~^ sec~^), MB trigger rate was typically 
10^° cm~^ sec~^ x 42 mb x 0.7 = 30 kHz. We have applied a prescale factor 
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to this trigger depending on its rate to keep the DAQ rate below an affordable 
level. Here a prescaled factor F is an integer number and defined as 

Number of triggers when the DAQ is alive 

F — — 1. 

Number of events recorded 

During the run, F = 19 to 79 were apphed to the MB trigger. 

2.7.2 Muon triggers 

Two kinds of muon-related triggers were prepared for physics triggers: 

(A) MB (8) ID for the single muon trigger and 

(B) MB IDIS for the dimuon trigger, 

where (8) stands for logical AND, "IDIS" and "ID" have been already ex- 
plained in section 2.5.3. 

The rate of the single muon trigger was usually less than 100 Hz and 
consistent with that of the irreducible hadron-decay background before the 
nosecone (1/1500 of the minimum bias rate) within a factor of two. How- 
ever, we sometimes observed abnormally high (~1 kHz) rates. This was 
found to be due to blown-up beams which hit the beam pipe and produced 
high-energy secondary particles (including muons) sailing through the MulD. 
(This background can be rejected with a vertex cut for offline analyses.) Wc 
used this trigger with F = 4 when the trigger rate was too high. Rates of 
the dimuon trigger were typically about 1/10 that of the single muon trigger 
and never beyond 100 Hz. Therefore no prescale factor was applied to the 
dimuon trigger that has been used for the J/ ip analysis described in the next 
section. 

Figure 16 shows Trigger Circuit (TC) efficiencies (etc) for 34 runs ran- 
domly selected from the runs used for the J/ip analysis. The efficiencies Etc 
reflect the inefficiency of the hardware trigger circuit. It is defined as 

number of events with both hardware and software triggers 
number of events with software triggers 

using MB-triggcred events, where a hardware trigger equals the NIM-logic 
trigger output (GLl input) and a software trigger is a result of the software 
which emulates the trigger algorithm. On average, an efficiency of 96.8 % 
has been obtained for the single muon trigger and 98.7 % for the dimuon 
trigger with small satistical errors (< 0.1%). The inefflciency for the dimuon 
trigger, 1.3%, is much smaller than other errors on the cross sections for 
production, which will be described in the next section. Inefficiencies 
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are ascribed to hardware dead time which depends on the trigger rate. A 
simple model calculation, in which all efficiency loss is assumed to be due to 
the dead time, reproduces Etc for the single muon trigger well, as shown in 
Fig. 17. This consistency ensures that both hardware and software triggers 
worked as expected. 



Trigger Circuit Efficiency 




IVIBxIDIS 
. — IVIBxID 



Figure 16: Trigger circuit efficiency (defined in the text) as a function of 
the run number for the single muon (MB ® ID) and dimuon (MB ® IDIS) 
triggers. 



2.8 Online Monitoring 

Conditions of the detectors in PHENIX have been continuously moni- 
tored online throughout the run to keep the quality of data. Monitoring 
of high-voltage (HV) status and data qualities are described in the follow- 
ing. In addition, status of such as low- voltage, electronics, and gas fiows was 
monitored, which was rather stable. 

High Voltage 

Status of all HV chains was monitored attentively, which were the most 
delicate component of a detector. Both the MuTr and MuID suffered from 
trip HV chains. They were sometimes due to high current caused by unstable 
beam condition, but usually due to spark current caused by humidity or other 
reasons. Flowing gas (air or nitrogen) into the secondary volumes of the 
detectors (just outside of the ionization gas volume) significantly reduced the 
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TC efficiency for MBx1 D 




I-- - Model 
■ — Measured 



Figure 17: Trigger circuit efficiency (defined in the text) as a function of tlie 
run number for tfie single muon trigger, witfi tfiat determined by a simple 
model calculation. 



number of trip chains for both the detectors. Trip chains were automatically 
recovered by a script monitoring and recording the status of all chains in 
every 10 (MuTr) or 60 (MuID) seconds. Because of shorter recovery time (1 
minute) compared to the typical run scale (1 hour) and small trip frequency 
1 trip per channel per run) for most of the channels, its effect on the 
variation of J/ip detection efficiency is minimal, which is conffimed in the 
next section. 

Data quality 

For both the MuTr and MuID, hit occupancies in each chamber plane 
were inspected for dead electronics, dead HV chains, and hot channels. For 
the MuTr, distributions of relative ADC counts of samples were monitored 
as shown in Fig. 18, to ensure the proper timing of the read-out (peaks in 
the third sample). For the MuID, TC efficiency, described in section 2.7, 
was checked to make sure the NIM-logic trigger was working properly. No 
significant deviation of these quantities from criteria was observed throughout 
the p+p run. 

2.9 Run and data summary 

Attached figure 17 shows the integrated luminosity RHIC has delivered 
to PHENIX in the Run-2 p+p period as a function of day starting from 
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Figure 18: Distributions of relative ADC counts, (ADC3 - ADC2) and 
(ADC3 — ADC4), which were monitored onhne to confirm their average 
values were positive. 

December 21, 2001 to January 23, 2002. It has continuously grown since the 
beginning of January and finally reached 700 nb~^. Due to finite live-time 
of the PHENIX DAQ system and the vertex cut of the BBLLl algorithm 
(75 cm), about a 150- nb~^ luminosity has been recorded with PHENIX. 
Table 5 summarizes the numbers of events with the triggers related to muon 
analyses. 



Trigger type 


Number of events 


Minimum bias 


196 xlO^ 


Single niuon 


31 xlO" 


Dimuon 


4.8 xlO^ 



Table 5: Numbers of events with the triggers related to muon analyses ob- 
tained in Run-2 p-|-p collisions. 

For the J/t/j analysis, the following run and event selections have been 
applied: (i) runs after RHIC has achieved stable beam condition and lumi- 
nosity (January 8, or the 19th day in Fig. ??), (ii) runs useful for physics 
analysis, excluding dedicated runs for beam and detector response studies as 
well as runs with HV or LV off, and (iii) events with a useful 2;- vertex, which 
is \zytx\ < 38 cm. 
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Finally, a 81-nb~^ luminosity has been used for the J/t/j analysis. The 
detail procedure to obtain it as well ets J/ip cross sections will be described 
in the next section. 



3 Data analysis 

In this section, the branching fraction for the decay J/^ — > {B^^) 
times rapidity-differential cross section for inclusive J/ if) production in the 
South Muon Arm acceptance B ^^da j / dy\y=i_-j is obtained, which is decom- 
posed into 



'MM 



dy 



J/ip particles are identified in an invariant mass spectrum of fi^fi^ pairs, 
where the number of J/ip^s, Nj/^, is counted with a reasonable background 

subtraction. The detection efficiency for the J/ip ^ events, s'H^ , is 

determined with a simulation, where detector response is well tuned to the 
real data. An integrated luminosity, £, is estimated also with a simulation, 
where consistency with the real data is confirmed. The rapidity width Ay is 
1.0 for this case. Details of analysis procedures to obtain A^j/^, and C 
arc described in the following. 

The rapidity range is further divided into two in each of which differential 
cross section is obtained (1.2 < y < 1.7 and 1.7 < y < 2.2) since distribution 
is expected to change significantly in this region. Similarly double-differential 
cross section in the invariant form, Bnfj^d'^aj/^/27rpTdpTdy\y=i,7, and the av- 
erage value of the transverse momenta of J/ip's, {pr), are obtained. 



3.1 Data sample 

The number of events with the triggers related to muon analyses have 
been shown in Table 5 in section 2.9. Dimuon triggered events, which contain 
most of the J/t/j yield (> 90%), have been used to obtain Nj/^. Minimum- 
bias triggered events and single-muon triggered events have been used to 
evaluate detector performance such as chamber efficiencies. 



3.2 Muon track reconstruction 

Muon tracks are reconstructed using offline software from raw data which 
contain information on MuID channel hits, MuTr cathode-strip charges and 
BBC ;2-vertices. This section is devoted to describe the muon reconstruc- 
tion algorithm, which starts from finding a track seed (a road) in the MuID 



38 



J/V' Production in p+p Collisions at ^/s = 200 GeV 



then grows by attaching clusters (consecutive cathode strips) at each MuTr 
station. 

3.2.1 MuID road finding 

First, one-dimensional roads are searched using either horizontal or verti- 
cal tubes. A road seed is constructed from a hit tube in the seed gap (or the 
first gap in the "search order" ) and event vertex position measured with the 
BBC or other interaction trigger counters. It is then extrapolated to the next 
gap in the search order, where additional hits are searched which are consis- 
tent with the road trajectory. The hit closest to the projection of the road to 
the gap within the search window is generally accepted and attached to the 
road. The road is fitted including the new hit and extrapolated to the next 
gap. The size of the search window is set to 15 cm (about two-tube widths) to 
allow a deviation from a straight-line trajectory due to multiple scattering in 
the steel absorbers. Two search orders, 2— >1— >3— >4— >5 and 3— >2^1^4— >5 
where each number n stands for the n-th gap of the MuID, arc used to 
enable roads to be reconstructed even the seed gap has inefficient tubes. 
Generally shallower gaps come earlier because of smaller multiple scattering, 
except for the first gap because of its higher hit occupancy from hadron and 
soft-electron background than the second or third gap. The algorithm is 
allowed to skip gaps without hits, in order to keep higher efficiency despite 
low-efficiency tubes in one or two gaps. Finally both horizontal and verti- 
cal roads are combined and final two-dimensional roads are reconstructed. 
Typical road-multiplicity is small (< 0.01) for minimum-bias events. 

Road finding efficiency is expected to be over 99% for muons with a 
momentum p > 3 GeV/c with good chamber efficiencies (97%) and low hit 
occupancies as in p-|-p collisions. However, it was actually lower than this 
since chamber efficiencies measured during Run-2 were lower (~90%), which 
will be described in the later section. 

3.2.2 MuTr cluster finding / fitting 

When a charged particle passes through a MuTr gap which is composed of 
two cathode planes and one anode- wire plane, charges are usually induced in 
two or three consecutive cathode-strips in each cathode plane. The hit posi- 
tion of the particle in each cathode plane is reconstructed with the following 
method. 

The amount of the peak charge in each cathode strip is determined by 
four ADC samples as described in section 2.5.1. A sequence of consecutive 
hit strips, called a cluster, is searched and fitted with an empirical formula 
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(Mathieson function [59]) to find the one-dimensional position at wfiicfi a 
particle would have passed in each cathode plane. Figure 19 shows an exam- 
ple of peak charge distribution for a cluster. The cluster position is obtained 
as the peak position of the fit function. 



S 120 
a> 

CI 

a 100 
.c 

o 

1i 80 

0) 
0. 

60 
40 
20 



-2-1012 

Strip Position (cm) 




Figure 19: An example of peak charge values of sequential cathode-strips 
(a cluster) induced by a charged track, together with a fit of the Mathieson 
function. In this example, the cluster has three hit strips and its position is 
determined to be 3.1-mm off the position of the central strip. 

If relative gain fluctuation and noise level arc 1% of a typical signal pulse 
or less, position resolution of 100 fim is obtained, which is confirmed with 
the cosmic ray test described in section 2.5.1. With this resolution, about 
110 MeV/c^ is expected as the mass resolution for a. J/ip. However the actual 
noise level has turned out to be worse during the run at station 2, which was 
typically 3%. Degradation in J/ip mass resolution due to higher noise level 
is expected to be about 30% (140 MeV/c^), which is compared with the real 
data in section 3.4. 

3.2.3 Muon track finding / fitting 

Starting from a road found in the MuID, a track grows by attaching 
clusters in the MuTr stations from backward to forward, that is, station 3 to 
station 1. 

Each road is extended to each cathode plane in station 3, where clusters 
are found consistent with the road inside its search window. Table 6 shows 
widths of the search windows for both the 9 and directions at each station in 
the South MuTr. The search window in station 3 is large enough compared 
to position resolution of roads (10 to 20 cm depending on quality cuts on 
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them). It gets smaller from station 3 to station 1 because track information 
such as position and momentum is getting more accurate in the process of 
track finding. By fitting clusters, a local vector called a stub is found. A stub 
requires at least two cathode planes with a cluster in each station. A crudely 
estimated momentum is assigned to the track from the last penetration gap 
of the MuID road. 



station 


6 direction (cm) 


(j) direction (cm) 


1 


20 


10 


2 


30 


25 


3 


50 


40 



Table 6: Widths of the search windows for the 9 and (f) directions at each 
MuTr station. 



The track is then extrapolated to station 2, using an effective bend-plane 
and a momentum kick, which are determined by the approximate momentum 
and the magnetic field inside the Muon Magnet. Clusters are searched again 
in station 2 within the search window. If a stub is found in station 2, the stubs 
in station 2 and 3 are fitted altogether to assign more accurate momentum 
to the track. It is again extrapolated to station 1 using better momentum 
information, where clusters are searched and attached to the track once more. 
The minimum number of hit planes required for a track is 12 out of all 16 
planes. 

Finally all the cluster positions and the event vertex position are fitted 
with the Kalman-Filter algorithm [60], which is a recursive technique to 

obtain the solution to a least-squares fit, to determine the momentum vector 
of the track at the event vertex taking energy loss in the absorber (the Central 
Magnet and copper nosecone) into account. 

Position resolution of the BBC, expected to be about 2 cm with a simu- 
lation, is confirmed with the real data. Figure 20 shows distribution of the 
differences between z-vertices found with the BBC and the Pad Chambers 
(PC) in the Central Arms [49]. The distribution is fitted with a double Gaus- 
sian function. The larger Gaussian has a 2.7 cm width, while the smaller one 
has a 9.8 cm width which is supposedly due to background contribution. 
Resolution of the PC is expected to be also about 2 cm with a simulation, 
which is consistent with the real data with background subtracted. 
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Figure 20: Distribution of the differences between 2;-vertices found with the 
BBC and the Pad Chambers (PC) in the Central Arms. The distribution is 
fitted with a double Gaussian function. The larger Gaussian has a 2.7 cm 
width, while the smaller one has a 9.8 cm width which is supposedly due to 
background contribution. 

3.3 Single muons 

In this section, some properties of reconstructed muons are shown to 
demonstrate their qualities and South Muon Arm performance. 

3.3.1 Event vertex distribution 

Low-momentum {p < 5 GeV/c) single muons, the majority of the inclu- 
sive muon yield, are expected to be dominated by charged hadrons decaying 
weakly into muons (tt^ — > /x^z/ and — > /J.^i') in flight before the nosecone. 
Contribution of punch-through hadrons is small, which has been confirmed 
by the beam test described in section 2.5.2. The decay probability for charged 
pions P(7r — > /i) is given by 

P(7r — >• /x) = 1 — exp(— Li?r/7CT) ~ LBr/'fcr 

where L is the distance between the vertex point and nosecone including 
one hadronic interaction length of the nosecone (15 cm), Br is the branching 
fraction for the decay tt^ f^"^^, 1 is the ratio of the particle energy to 
mass and cr is the decay length of the particle. Typical value of LB^/^ct 



42 



J/V' Production in p+p Collisions at ^/s = 200 GeV 



is 3xl0~^ for p ~ 3 GeV/c muons. The same kind of formula holds also for 
charged kaons. The net contribution of charged kaons to single muon spectra 
is expected to be about the same as that of charged pions. 

Figure 21 shows BBC z-vertex (z^utx) distribution for single-muon events 
divided by that for minimum-bias events in the range \zytx\ < 38 cm which is 
a cut used for the J/i/j analysis. Error bars indicate statistical errors of the 
particular run. The solid line shows a straight line fit to the data assuming 
flat distribution for non-decay components, which reproduces the data very 
well in the range —20 < z^tx < 38 cm. The deviation from the fit in the 
range —38 < z^tx < —20 cm is due to the increase in background induced by 
very low angle particles. From the fit, the fraction of non-decay components 
is determined to be about 20%, which supposedly includes background such 
as punch-through hadrons and ghost tracks as well as physics signals such as 
charm and bottom mesons decaying semi-leptonically. 

The dominance of decay muons to the measured single-muon yield is 
consistent with an expectation, which confirms that the South Muon Arm 
has been worked as expected and there is no significant contribution of ghost 
tracks nor punch-through hadron background to reconstructed muons. 




Figure 21: BBC 2;- vertex distribution for single-muon events divided by that 
for minimum bias events. Error bars include statistical errors only. The solid 
line shows a linear fit to the data points assuming constant distribution for 
non-decay components. 
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3.3.2 Track-road matching 




Figure 22: The distribution of the distance between a track and a road 
(T^track-road) ^t station 3. The peak position, around 12 cm, roughly represents 
the rms in the two-dimensional space. 

Figure 22 shows distribution of the distance between intersections to a 
station-3 plane of a track and a corresponding road, rtrack-road- The peak 
is around 10 cm which is consistent with the expectation from the position 
resolution of a road (8.4 cm) and multiple scattering in the absorbers, thus 
demonstrating that the two detectors match well as expected. 

3.4 Dimuon mass and A^j/^ 

Figure 23 shows invariant mass Minv spectra for both unlike-sign and 
like-sign muon pairs with the following cuts: 

• BBC vertex position |2;i,t-r(BBC) | < 38 cm, which is between the north 
and south copper nosecones (|2;| < 40 cm) minus resolution of the BBC 
(2 cm), and 

• track vertex consistency \zvtx(!^^)—Zvtxi!^'^)\ < 30 cm where ^^^(Trl) 
and ^„ta;(Tr2) represent the vertex position of each track obtained as 
the closest approach to the 2;-axis. This cut corresponds to about 2a 
of the position smearing of a track due to multiple scattering in the 
Central Magnet steel (about 10 cm for each track). 
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Figure 23: Invariant mass spectra for unlike-sign (points with statistical 
errors) and like-sign (line) muon pairs. Cuts are described in the text. 



A significant enhancement is found in the J/ip mass region for unlike-sign 
pairs while no such a peak is found for like-sign pairs. 

Figure 24 and 25 show BBC z-vertex distributions for unlike-sign muon 
pairs in the J/t/j mass region (2.5 GeV/c^ < Minv < 3.7 GeV/c^) and the 
background region {Minv < 2.5 GeV/c^) respectively. Larger contribution 
from the far-side of the South Arm (positive z) in the background region 
confirms that the dominant fraction of background comes from hadron decays 
as shown also for inclusive single muons (Fig. 21). On the other hand, the 
flat distribution for the J/ip candidates is consistent with the much shorter 
decay-length of a J/-^ (cr = 2.3 x 10~^^ m). 

From a simulation study using the PYTHIA event generator [61] and 
GEANT [58], it is found that the combinatoric hadron-decay background 
docs not produce a large (~30%) difference between the numbers of unlike- 
sign and like-sign muon pairs in the J/ip mass region including dead elec- 
tronics and chamber efficiencies during the run. The consistency of PYTHIA 
simulation with experimental data will be examined in section 3.6.3. The 
number of fi^fi^ pairs from the Drcll-Yan process is estimated to be small 
(~0.03 events) according to a PYTHIA simulation. Thus we have simply 
assumed that the number of background for unlike-sign pairs is the same as 
the number of like-sign pairs then obtain 

Nj/^ = 41 (unlike-sign pairs) — 5 (like-sign pairs) = 36 ± 6.8(stat.). 

Figure 26 shows the difference between the invariant mass spectrum of 
unlike-sign muon pairs and that of like-sign pairs, which is fit to a Gaussian 
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Figure 24: BBC z-vertex distri- 
bution for unlike-sign muon pairs 
with 2.5 < Mir,o < 3.7 GeV/c^. 
The hne shows a straight hne fit 
to the data. 
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Figure 25: BBC z-vertex distri- 
bution for unhke-sign muon pairs 
with Minv < 2.5 GeV/c^. The 
line shows a straight line fit to the 
data. 
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Figure 26: The difference between the invariant mass spectrum of unlike-sign 
muon pairs and that of like-sign pairs, which is fit to a Gaussian distribution. 
Cuts are described in the text. 
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distribution. The center position, 3156 ± 74 MeV/c^, is consistent with the 
J/ijj mass (3097 McV/c^ [62]). The width, 257 ± 75 McV/c^, is shghtly worse 
than the expectation with the reahstic chamber noise (3% at station 2), which 
is roughly 140 McV/c^. This discrepancy is supposedly due to the alignment 
of MuTr chamber positions done only to 300 /xm due to the limited number 
of straight tracks with the magnetic field off used for the alignment. Mis- 
alignment effectively worsens the chamber resolution since it can differ octant 
by octant and has different effects on each charge. The net effect on the 
J/ip mass resolution is studied with a simulation and about 190 MeV/c^ is 
obtained as shown in Fig. 27, which is consistent with the real data. 



J/\\f mass (simulation) 



h101 

Nent = 299 
Mean = 3.093 
RMS =0.1943 

Chi2/ndf = 13.99/17 
Constant = 35.07± 2.702 
Mean = 3.093± 0.01138 
Sigma = 0.1866+ 0.00902' 




Dimuon invariant mass in GeV/c 



Figure 27: Invariant mass spectrum for unlike-sign muon pairs from J/?/) — > 
jJ.'^IJi'' events obtained with a simulation. MuTr cathode planes are mis- 
aligned intentionally by about 300 fim.. The curve shows a fit with the 
Gaussian function. 



We estimate a systematic error of 10% on the background-subtracted yield 
adjusting the mass window by 300 MeV/c^ that is roughly the resolution for 
a J/V'. The Gaussian fit gives Nj/^ — 33.0 which is consistent with the 
simple counting method within the systematic error. 

To obtain pt and rapidity differential cross sections, invariant mass dis- 
tributions for two rapidity bins {1.2 < y < 1.7 and 1.7 < y < 2.2) and four px 
bins (0 < < 1, 1 < Pt < 2, 2 < < 3 and 3 < pt < 5 in unit of GeV/c) 
are shown in Figs. 28 through 29 and Figs. 30 through 33 respectively. The 
numbers of J/ip^s in each bin are obtained in the same method described 
above and summarized in Tables 9 and 10 in section 4.1. Fluctuations of 
the background yields in each range seem to be of statistical nature and no 
effect is expected on Nj/^. For the highest px bin (3 < pr < 5 GeV/c) with 
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Nj/ip — 3, background contamination due to such as mis-alignment of MuTr 
cathode strips is estimated to be small (<^1 count) from the counts in the 
higher mass and higher pt regions as summarized in Table 7. 





2.5< Mj/^ <3.7 GeV/c2 


Mj/^ >3.7 GeV/c2 


3< pt < 5 GeV/c 


3/0 


0/1 


Pt > 5 GeV/c 


0/0 


2/2 



Table 7: Dimuon counts in high-p^ and high mass regions. Numbers of 
both unlike-sign pairs (left-hand side) and like-sign pairs (right-hand side) 
are shown in each column. The underline represents that the number is for 
the signal. 
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Figure 28: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs in 1.2 < y < 1.7. 



Figure 29: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs in 1.7 < y < 2.2. 



3.5 J/i/j detection efficiency 

Total detection efficiency for J/ijj ^ A*"^/^" events in p-|-p collisions, s///, 
is decomposed into four factors: 

J/'>P _ . ■l/^'P . J/'>P . J/'>P (r,\ 

^tot — 'lace ^MuID ^MuTr ^ BBC \^ ) 

where 
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Figure 30: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs with < pr < 1 GeV/c. 



Figure 31: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs with 1 < < 2 GeV/c. 
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Figure 32: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs with 2 < < 3 GeV/c. 



Figure 33: Invariant mass spectra 
for unlike-sign and like-sign muon 
pairs with 3 < pt < 5 GeV/c. 
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South Muon Arm acceptance times reconstruction efficiency 
for muon pairs from J/ip^s produced in 1.2 < y < 2.2 with 
100% chamber efficiencies, 

Efficiency correction due to real chamber efficiencies of the 
MuID, 

Efficiency correction due to real chamber efficiencies of the 
MuTr, and 

Efficiency of BBC for p+p ^ J/^jX (1.2 < y^/^ < 2.2) events. 

Analysis procedures and results of each factor will be described in the fol- 
lowing. 

3.5.1 Detector acceptance 

The geometrical acceptance of the South Muon Arm is represented as rjacc 
including reconstruction efficiency for muon pairs from J/it''s produced in 
1.2 < y < 2.2 with 100% chamber efficiencies. In p+p events, multiplicity is 
sufficiently small to achieve high (~90%) reconstruction efficiency when hits 
are found in chambers. Therefore, rjacc is close to the geometrical acceptance 
itself. Using simulation, rjacc is calculated as 

the number of reconstructed J/ijj events with 1.2 < y < 2.2 

rjacc ~~ 

"'^'^ the number of simulated J/ip events with 1.2 < y < 2.2 

where PYTHIA has been used to produce J/ip events, based on the color- 
singlet model. Rapidity and pt distributions (at lower px) of J/ip do not 
depend on the production model unlike a total cross section, and are consis- 
tent with the real data as discussed in section 4.1. The same reconstruction 
software, parameters and cuts are used as the real data. The result of the 
average value of rjacc is 0.11 requiring a dimuon trigger. Rapidity and pr 
dependence will be described in section 3.5.5 including finite chamber and 
BBC efficiencies (2^/). 

Polarization of J/i/j or spin- alignment (denoted as A), which is sensitive 
to its production mechanism, is unknown at ^/s = 200 GeV and not possi- 
ble to be determined with the limited number of J/ijj^s obtained in Run-2. 
The South Muon Arm acceptance has significant A dependence as shown 
in Fig. 34, giving a systematic uncertainty of rjacc- This dependence is pri- 
marily a result of the pz < 2GeV/c lower-momentum cutoff which cuts off 
the backward-going daughter muons in the J/ip rest-frame with respect to 
the momentum direction of the J/ijj in the laboratory frame. Results of both 
lower-energy experiments and Tevatron indicate that |A| is no larger than 0.3 
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especially for low-p^ (~ 1 GeV/c) and low-xp (~ 0.1) J/ip^s [43,47] which 
dominate our yield. On the assumption of |A| < 0.3, we have assigned a 10% 
systematic error of rjacc- 
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Figure 34: J/ip polarization dependence of the J/i/j acceptance at pt = 1 
GeV/c and some y (rapidity) points. 




3.5.2 MuID efficiencies 

MuID efficiencies for single- muons and dimuons are estimated using MulD- 
panel efficiencies obtained from a sample of 37 runs randomly picked up to 
cover the entire run-period. 

Efficiency of the specific MuID plane d for a muon, Ed, is given by 

the number of reconstructed roads with a hit in d 
the number of all reconstructed roads 

for each orientation. Only minimum-bias triggered events have been used for 
this efficiency study since MuID-triggered (single-muon and dimuon trigger) 
events are possibly biased by the trigger algorithms. Since we allow "skipped 
gaps" in the road finding algorithm described in section 3.2.1, the plane 
d is not necessarily excluded from the search order. This is confirmed by 
comparing results with two search orders with and without the plane d. Since 
efficiencies for planes in the seed gaps^ are not possible to be extracted, two 
sets of search orders are used to obtain panel efficiencies for all the gaps: 



The first and second search-order gaps in this case. 



R D. Sato 



51 



• search order set A: 2^1— >3— >4— >5 and 1^2— >3— >4— >5 for gap 3,4,5 
and 

• search order set B: 4— >5— >3— >2— >1 and 5— >4— >3— >2— >1 for gap 1,2,3 

where two search orders in each set are tried in the algorithm. Panel efficien- 
cies in gap 3 obtained with A and B are compared and found to be consistent 
within their statistical uncertainties (2%), which indicates panel efficiencies 
are not sensitive to the search order. 

Tighter cuts in the road finding algorithm, such as the rz=o (the distance 
between the origin and the intersection of a road to the z = plane) cut, are 
used than those for the usual muon reconstruction to reject ghost roads as 
much as possible and select only good quality roads. Furthermore, additional 
cuts are applied to sample roads: (i) each road should penetrate through all 
the five gaps, (ii) only one road is found in each event, (iii) each has a hit in 
the other orientation in the panel of interest and (iv) each has a hit of either 
orientation in both adjacent planes. Cut dependence on panel efficiencies is 
studied for example, change (iv) (iv)': has a hit of both orientations in 
both adjacent planes. The result is consistent with the original one to the 
statistical accuracy of the result (2%). 

Figure 35 shows distribution of panel efficiencies for all the 60 panels in the 
South MuID including both orientations. Most of them gather around 90% 
but some are smaller because HV (high-voltage) potentials applied to both 
layers of each panel-orientation (see Fig. 13 in section 2.5.2) were not high 
enough (< 4300 V) in order to avoid sparks or some problems. In addition, 
about one third of all HV chains had one dead or lower-potential layer, for 
which about 90% efficiency is expected. However, other two thirds had good 
potentials on both layers, where about 98% efficiency is expected [57]. The 
panel efficiencies obtained from the real data are somewhat smaller than the 
naive expectation. A possible explanation is we have applied potentials of 
4350 V nominally during the run which may not be high enough since it is 
about the edge of the plateau of the excitation curve (typically, 4300 V). For 
example, an effective potential may drop by 100 V with a 100 nA current 
fiow on an anode wire, where a typical measured-current on one HV chain 
(about 20 tubes) was 10 to 100 nA. 

Trigger and road-finding efficiencies have been estimated with simulation 
including the trigger emulator and the MuID chamber efficiencies obtained 
above for both single muon and dimuon (J/i/j) events. Definition of MuID 
efficiency for single muons SmuID given by 
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Figure 35: Statistics of MuID panel efficiencies for a single muon. 



where represents the number of reconstructed muon tracks through a full 
simulation. The similar definition is applicable for J/V^'s, which is equivalent 
to £muid requiring a dimuon trigger in both the numerator and denomina- 
tor. For a single muon with Pz — 5 GeV/c, 80% efficiency is obtained and 
62% for a J/'^ requiring a single or dimuon trigger respectively. These are 
consistent with a very crude prediction that e'l^^j^, ~ {^^muId)'^ since the av- 
erage momentum from J/'^'s (p^ = 5 GeV/c) is used for the single muon 
simulation. 

Without trigger requirements higher efficiency is expected, since the of- 
fline road-finder is more flexible than the online trigger algorithm, that is, 
less susceptible to chamber efficiencies. For single muons, 86% efficiency is 
obtained using minimum-bias triggered events. The ratio {£muID requiring 
a trigger) to {s^m^jd requiring a trigger) is obtained, which is 0.93 and 
consistent with that of the real data (0.92) with about 2% statistical errors 
for each. This consistency verifies the chamber efficiencies obtained with the 
method described above. 

Run dependence of £^^7^ is also studied with a simulation and found 
that its variation is 11% which is shown in Fig. 36. This fluctuation is 
roughly consistent with statistical uncertainties of panel efficiencies which are 
typically 3% for each panel and the minimum requirement for the number of 
hit planes which is 9 for a dimuon trigger (3% X\/9 = 9%). 
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Figure 36: Run by run fluctuation of MuID efficiencies for J/ijj^s (^mu/d) 
Dimuon triggers are required. 



MuID Hit Occupancy 



In high particle-multiplicity environment such as in Au+Au collisions, 
road finding efficiency may drop because of difficulty in hit-road associations. 
However in p+p collisions, much smaller hit occupancy ensures high road- 
finding efficiency. Figure 37 shows distribution of the number of MuID hit 
channels per event. Most of the minimum-bias triggered events have zero 
or small hit multiplicities. On the contrary, the hit multiplicity distribution 
for the single- muon triggered events has a peak around 9 or 10 which is the 
expected number of hits when a single muon penetrates through the entire 
MuID (5 gaps X 2 orientations). This indicates that most of MuID hits 
are induced by signal muons and not by background hits such as from soft 
electrons. 

Figures 38 and 39 also show MuID hit multiplicities in single-road events 
and J/ip candidate events (unlike-sign dimuons in the J/ip mass region) 
respectively, both of which are compared with simulations. For both cases, 
hit multiplicities agree with simulations and again confirm hit occupancies 
in the real data are small as expected. Provided valid particle multiplicity 
of the event generator (PYTHIA) which is confirmed in the next section, 
this consistency assures response and performance of the detectors are as 
expected; thus, no efficiency degradation is suspected. 
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Figure 37: MuID hit multiplicity in minimum-bias triggered events (dotted 
line) and single-muon triggered events (solid line). 




Figure 38: MulD hit multiplicity for single- road events. The solid line shows 
real data and the dotted line shows simulation data. 
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Figure 39: MuID hit multiplicity for J/ip (candidate) events. The sohd hne 
shows real data and the dotted line shows simulation data. 

3.5.3 MuTr efficiencies 

Efficiencies of each MuTr plane for a muon are known to be high (99%) 

inside active detector volumes according to a cosmic ray test. Inactive vol- 
umes between octants, which are their frames, are included in the detector 
acceptance (about 70% for a single track). Some dead Front-End Modules 
(FEMs) and high-voltage supply chains during the run caused inefficiencies 
primarily. 

A sample run was selected to model MuTr dead FEMs, dead electronics 
channels and dead high- voltage chains. To examine consistency of the model 
with the real data, MuTr efficiencies for a single muon defined as 

„ the number of roads with a track 

the number of roads 

are obtained using both the real data and simulation to be compared. Tighter 
cuts are applied to the sample roads to reject as many ghosts as possible. 
Figure 40 shows e^muTt ^ function of r^=0; which is the distance between 
the origin and the intersection of a road to the z = Q plane. The decrease 
in efficiency at large rz=Q is due to contamination of background roads pre- 
sumably from beam scraping that has caused high trigger rates as described 
in section 2.7.2. The cut rz=Q < 50 cm is applied to eliminate them. Other 
road quality parameters are checked in the same way and confirmed to be 
valid. About 80% of the road sample is dominated by hadron-decays as 
demonstrated by the BBC vertex distribution as already shown in Fig. 21. 
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A fraction of ghost roads is further constrained by z-vertex dependence of 
^'muTti because efficiency is sensitive to the fraction of non-ghost roads with 
2;-vertex dependence, which is dominanted by decay muons. Actually, there is 
no significant dependence on BBC 2;-vertex as shown in Fig. 41. The fraction 
of the ghost roads to the road sample is constrained by fitting it assuming a 
flat z-vertex distribution for them and estimated to be less than 5%. 
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Figure 40: MuTr efficiency SmuTt ^ function of r^=o of roads. Sec the 
text for the definitions of ^MuTr ^z=o- The dotted line and the arrow 
represent the quality cut to select good roads used for the MuTr efficiency 
calculation. 



Figure 42 shows MuTr efficiencies for both the real and simulation data as 
a function of the azimuthal angle of the qualified roads. There is a reasonably 
good agreement of the real data and simulation with exception of the region 
in = —45° to —22.5° (octant 8, half-octant 1 or 8-1). This discrepancy is 
turned out to be due to the improper treatment of the simulation data in 
8-1, which has reduced its acceptance significantly (25%). Since an opening 
angle between two muons from a J/ip decay is large (30 ~ 60 degrees), the 
probability for either muon to fall into the half-octant 8-1 is roughly (1/16) 
X 2 = 1/8 where 16 is the number of all half-octants. In that case, efficiency 
is 25% of the ordinary case. Therefore, the overall efficiency for a J/i/j, ^'muTt 
is estimated to be 

(^1 - + ^ X 0.25 ~ 0.90 

compared to the normal case. Therefore we assign a systematic error of 10% 
to £MuTr ^bis particular run based on the disagreement between the sim- 
ulation and real data. The 0-averaged values of EmuTt obtained are (29.6 
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Figure 41: MuTr efficiency ^^muTt ^ function of the BBC z-vertex. See the 
text for the definition of s^muTv "^^^ error bars are statistical errors for one 
particular run. The line shows a ffi to the data with a function assuming a 
constant fraction of ghost roads, which is determined to be less than 5%. 




Fi gure 42: MuTr efficiency s^^^r ^ function of the azimuthal angle of 
roads. See the text for the definition of s^i^^xr- '^^^ points with a statistical 
error are for the real data and the histogram is for the simulation data. 
Statistical errors of the simulation data are omitted, which are comparable 
with those of the real data. 
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± 0.3) % for the real data and (28.2 ± 0.3) % for simulation with statisti- 
cal errors, whose difference is consistent with the systematic error of SMuTr 
obtained above assuming a crude relation s'^^^j,^ ~ {^MuTrY- 

To minimize an uncertainty due to run dependence of MuTr efficiency, 
only runs with a high duty fraction are selected and used for the analysis. 
Here we define the duty fraction DF for each run as 



the number of chains with high- voltage on 
the number of all chains 



each integrated over all measurements in every 10 seconds during the run. 
We have selected runs with DF > 0.89. To estimate run dependence of MuTr 
efficiencies, 26 sample runs with DF > 0.89 and relatively longer run time 
are selected to see their variations. Figure 43 shows distribution of s'^^^j.^ for 
those runs and variation is found to be 4%. It should be noted that e^uTr 
slightly higher than (^^/ut^^)^ ~ 9% since the average momentum of muons 
from J/ip decays (about 5 GeV/c) is higher than that of muons used for the 
calculation of SmuTt (about 3 GeV/c) with lower efficiency. 
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Figure 43: Run statistics of MuTr efficiency for J/ijj^s {£MuTr) f*^^ ^6 sample 
runs. 



Cluster hit occupancies in the real data have been examined. Figure 44 
shows distribution of the average number of clusters in each half-octant for 
minimum bias events. They are typically 0.02 or 0.03 and sufficiently small 
not to worsen reconstruction efficiency for muon tracks. 
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Figure 44: Distribution of the average numbers of MuTr clusters in each half 
octant for minimum-bias events. 



3.5.4 BBC efficiency for p+p J/ipX events 

Efficiency of the BBC for p+p — > J/ipX events {s'^Jbc)^ where a, J/ijj 
decays into a muon pair detected in the South Muon Arm, is estimated to 
be 0.74 ± 0.01 (stat.) from simulation studies using PYTHIA and GEANT. 
No significant pt nor rapidity dependence has been found. This factor is 
closely related to £^sc which is BBC efficiency for p+p inelastic events and 
will be described in section 3.6.3. The systematic error of £^bc '^^^^ ^® 
discussed there. 

3.5.5 Total detection efficiency 

Total detection efficiency e'Igf is given by 

j/^ the number of reconstructed J/^'s in 1.2 < y < 2.2 
*°* the number of simulated J/V^'s in 1.2 < y < 2.2 

with a simulation including the MulD, MuTr and BBC efficiencies obtained 
in the previous subsections. Run-averaged values are used for the MuID 
part. For the MuTr part, dead FEMs, high-voltage chains and real ADC 
gains with the particular run resulting 15.3% efficiency (see Fig. 43) are 
used. About 88 000 PYTHIA events with unpolarized J/ijj^s are simulated 
and e/// is obtained which is (1.19 ± 0.03)%. where the error stands for the 
statistical error of the simulation. PYTHIA's pt and rapidity distributions 
are consistent with the real data as shown in the next section. 
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Transverse momentum (pt) and rapidity (y) dependence are shown in 
Figs. 45 and 46 respectively. Error bars represent statistical errors of the 
simulation. Relatively small pt dependence of 5/// is found in the range 
< Pt < 5 GeV/c which is the statistical limit of our measurement. For the 
calculation of differential cross sections for J/ip production, averaged values 
are used in each rapidity or pt bin shown in Tables 9 and 10 in section 4.1. 
A slight decrease in s^^^ at higher pt is due to the bias of the dimuon trigger 
which requires two different quadrants to be fired, thus losing some J/ip^s 
with a smaller opening angle between two muons. A systematic error of e^lf 
due to the uncertainty of J/ip^s pt distribution is estimated to be small (< 
2%) by changing efficiency values in each px bin within their statistical un- 
certainties. Rapidity distribution is sensitive to gluon density in the proton. 
A systematic error of e/// due to the uncertainty of the rapidity distribution 
is studied using PYTHIA with various parton distribution functions. As a 
result the variation of e/// is found to be negligible (< 2%). 
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Figure 45: Transverse momentum (pt) dependence of J/ip detection effi- 
ciency {s'tot')- Error bars are for statistical errors of the simulation. 

The systematic errors of each factor of e/// is summarized in Table 1 1 in 
section 4.1. 



3.6 Integrated luminosity 

The integrated luminosity used for this analysis is given by 
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Figure 46: Rapidity dependence of J/ijj detection efficiency {sUf)- Error 



bars are for statistical errors of tfie simulation. 



^ B BC '^inela 

where 



Ar|2i,ta;|<38cm 

C = (3) 



^I.„^.l<d«cm. number of minimum bias triererers with an offline 

BBC 2;- vertex cut, 
Srbc" efficiency of the BBC for p+p inelastic events with the 

same vertex cut, and 
Cineia' P+P inclastic cross section. 



Analysis procedures to obtain these factors will be described in the following. 



3.6.1 Number of minimum bias triggers 

The number of minimum-bias triggered events integrated over all the 
selected runs is represented as N^^^^^^^^'^™, with an offline BBC 2;-vertex cut 
l^vtxl < 38 cm. Since most of minimum bias triggers have been prescaled (see 
section 2.7.1), it is not possible to obtain un- prescaled values of 
Instead it is calculated as 

nrrurii , | | <38 cm 

y.T-r-«nf,|2;„tx|<38cm _ Arruni MB, prescaled / ,\ 

-''MB ~ -''MB, live ^ ivrrurii 

-''MB, prescaled 
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for each run where un-prescaled (hve) counts of all minimum-bias triggers 
-^Ms'iive obtained from the run control log and the ratios of minimum bias 

events with a BBC ^-vertex cut (^I/i^' prcsciiecf '^™) minimum bias events 
(^MB prcscaicd) ^'^^ each run are obtained from the prescaled data. Intergated 
over all the runs used for the J/ijj analysis, 1.72 xlO^ has been obtained with 
a very small statistical error. 

The fraction of the number of triggers from beam-related background is 
found to be less than 0.1% with the beam test, where the rate of the BBLLl 
trigger was reduced from 10 kHz to 2~3 Hz when beams were mis-steered. 
Crossing-by-crossing variation for minimum-bias triggered events, shown in 
Fig. 47, also give us an estimation of the magnitude of the background. Since 
beam-related background events happen randomly, its fraction is estimated 
by the ratio of the number of events in non-collision crossings (even numbers 
in the figure) to the number of events in collision crossings (odd numbers in 
the figure). If a BBC 2;- vertex (|-2,,ja:| < 75 cm) is required (the upper figure), 
the ratio is less than 0.1%. On the other hand, it goes up to 10% if no BBC 
2;- vertex required (the lower figure), which means NTC-exclusive triggered 
events suffer from beam-related background by 10%. 

3.6.2 p+p inelastic cross section 

The total p-|-p cross section atot at ^/s — 200 GeV has not been measured 
yet. Following parametrization for <Ttot as a function of ^/s is given in [62] as 
well as fit parameters, which reproduce experimental data well. 

atot ^Z + B \og\s/so) + Y,{s)-^' - ^2(5)"^^ (5) 

where Z, B, are in mb and s and Sq are in GeV^. The first and second 
terms, which represent the pomeron contribution, dominates at high y^. 

The exponents r/i and 772 represent lower-lying C-even and C-odd exchanges, 
respectively. From the fitted parameters given there, atot is determined to be 
51.5 ± 1.1 mb at ^/s = 200 GeV. Elastic cross section aeia is also obtained 
(10.2 ± 0.3 mb), yielding the inelastic cross section of (Tj„e/a — (^tot ~ CTeia — 
41.3 ± 1.2 mb. 

Pseudo-rapidity range of protons from elastic p+p scattering is 5 < l^^l < 
9 which is out of BBC acceptance, thus being excluded from the luminosity 
calculation. 

3.6.3 BBC efficiencies 

Efficiencies of the BBC for inclusive p+p inelastic events, e^bbci 
p+p J/il)X events, Sbbc-> sensitive to charged particle multiplicity 
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Bunch-Crossing Distribution with BBC z-vertices 
















20 40 60 


80 100 



Bunch Crossing Number 
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Figure 47: Crossing-by- crossing variation for minimum-bias triggered events 
with (upper) and without (lower) a BBC vertex cut {\zytx\ < 75 cm). Even 
numbers are for non-coUision crossings and odd numbers are for collision 
crossings. Events in non-coUision crossings are supposedly induced by beam- 
related background. 
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which has not been measured at ^/s — 200 GeV before. We have used the 
event generator PYTHIA to estimate them, whose consistency with the real 
data is fully examined. 

PYTHIA and the PHENIX detector simulation using GEANT give e-^J^^ 
= 0.74 ± 0.01 and e'^^^^ = 0.51 ± 0.01 respectively with a \zytx\ < 38 cm 
cut, where errors are statistical errors of the simulation. Because of higher 
multiplicities in p+p — > J/ipX events than in inelastic events on average, 
^BBC larger than £^^c- significant pt nor rapidity dependence has 
been found for e'^^Q. 

To cross-check these results for their validities, relative efficiency R = 
^'bbc/^^mb^ compared with the real data, where is efficiency of the 

MB (minimum bias) trigger, logical OR of BBG and NTG triggers, for p+p 
inelastic events. Prom simulation, R — 0.51/0.70 = 0.73 is obtained with a 1% 
statistical error. However, NTG trigger rate has suffered from beam-related 
background which is not included in the simulation. Its fraction is about 
10% at maximum as described in section 3.6.1. Therefore R is expected 
to vary from 0.67 to 0.73. Figure 50 shows run statistics for the relative 
rate of the BBG trigger to the minimum-bias trigger which is identical to 
R. The distribution is consistent with the expectation from the simulation 
and background rate. The lower tail is due to nms with higher background 
(or NTG trigger) rate. Even at the highest rate, contamination of the beam 
scraping background to physics events is expected to be small (< 10~^). 

BBG efficiencies are sensitive to particle multiplicities especially in and 
near the BBG acceptance (3.0 < \r)\ < 3.9) given by the event generator. 
Pigure 48 shows comparison of PYTHIA with the UAl data [63] for the pt 
spectrum for inclusive charged hadrons at central rapidity {\ri\ < 2.5) in p+p 
collisions at y/s = 200 GeV. They are in good agreement to 20% precision. 
Since gluon-gluon scattering dominates the inclusive charged-particle yield 
for both p+p and p+p collisions, the consistency in p+p collisions indicate 
the consistency in p+p collisions. In Run-2 at RHIG, Pt differential cross 
sections for the neutral pion production have been actually measured in the 
PHENIX Gentral Arms (|?7| < 0.35) in p+p collisions at ^/s = 200 GeV which 
agree with the UAl results in the range < 4 GeV/c to about 10% precision 
[64]. Pigure 49 shows comparison of PYTHIA with the UA5 data [65] for 
the charged-particle pseudo-rapidity distributions in p+p collisions at ^/s = 
200 GeV. PYTHIA distribution is again consistent with the real data to 20% 
precision. A systematic error of e'^^^q from uncertainty of particle distribu- 
tion is estimated by changing charged particle multiphcity in PYTHIA by 
20% tuning parameters on the fragmentation function ^ and changes of 



5for example, change PARJ(41) from 0.3 to 1.3 and PARJ(42) from 0.58 to 1.58, which 
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are found to be 15% or less. PDF (parton distribution function) dependence 
is also studied and is found to be smaller (3%). Therefore a 15% system- 
atic error is assigned to s'bbc from the uncertainty of the initial particle 
distribution. The relative trigger rate R is 0.80 when particle multiplicity is 
increased by 20% and 0.66 when decreased, both of which are deviated from 
the measurement and indicating the particle multiplicity with the original 
PYTHIA parameters reproduces the real data well. 




Figure 48: Inclusive cross section for single charged hadrons in p+p collisions 
at ^/s = 200 GeV as a function of Pt of hadrons. The points with error bars 
are for the UAl data [63] and the histogram is for a PYTHIA prediction. 

In the p+p run, event vertex distribution varied from 40 cm to 70 cm 
in rms. If £^5^ had a significant dependence on the ^-vertex position, it 
would be affected by the change in the vertex distribution. However, no z- 
vertex dependence within \z^tx\ < 38 cm is actually found over the statistical 
uncertainty of the simulation study (< 3%); hence, no fluctuation in s^^^q is 
expected during the run. 

High-voltage potentials apphed to each PMT (photomultipher tube) of 
the BBC were stable during the run. However there is PMT by PMT fluc- 
tuation in discriminator-threshold values possibly by 10 to 20% which is not 
reflected in the simulation. It gives rise to an uncertainty of s'^^qq. However 
it is found to be small (< 5%) with a simulation by conservatively doubling 
and halving all the PMT threshold values. 



represent the parameters a and b respectively in the Lund fragmentation function f{x) oc 
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Figure 49: Charged particle pseudo-rapidity distribution in p+p collisions 
at ^/s — 200 GeV. The points with error bars are UA5 data [65] and the 
histogram is for PYTHIA. 



As the ver der Meer scan ^ result, e^^'c" ^ (^ineia has also been determined to 
be 18.5 ± 1.9 mb. This is consistent with e'^^^^ obtained with the simulation 
times aineia obtained in the previous subsection within their errors. 

The same magnitude of uncertainty is expected on e'^^^j. However, some 
systematic uncertainties of s^^^^, and e'^^^ will be canceled out since the 
total cross section is proportional to their ratio. For example, if particle 
multiplicity is higher in inelastic events, it tends to have higher multiplicity 
also in J/-^ — > n~^fi~ events. Therefore a smaller systematic error of s'^^q is 
assigned which is 10%. 

BBC efficiency for another hard scattering events in p+p collisions, 7r° 
production events, has been measured in the Central Arms (|7^| < 0.35) in 
Run-2 [64]. It is relatively px independent and consistent with a PYTHIA 
simulation within a few percent of errors, as shown in attached figure 18. It 
confirms that PYTHIA reproduces particle multiplicities in p+p hard scat- 
tering events for which BBC efficiency is scale and (supposedly) process in- 
dependent. 



technique to measure the cross section of a beam by steering 
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Figure 50: Run statistics for the relative trigger rate, BBC/MB, where MB 
(minimum bias) is the logical OR of the BBC and NTC triggers. The lower 
tail is due to the runs with higher NTC trigger rate from background. 



3.7 Average px 

The average value of the transverse momentum (p^) of J/'0's, (pt) is 
calculated as 

(Pt) = CeffCrnom-offset{PT) uncorr 

where {pT)uncorr IS an uncorrected value of (pt), Ceff is a correction factor 
due to Pt dependence of efficiency and C mom-offset is another correction 
factor due to the shift and resolution of momentum. Following descriptions 
are for procedures to obtain each factor. 



3.7.1 Uncorrected (pr) 

All data points arc < 5 GcV/c, which is reasonable with the num- 
ber of events obtained. A statistical error of (pr), ^{Pt) is calculated as 
V {Pt) /-^ 1 where N is the number of events. Table 8 shows (pt) and A(pr) 
for both unlike-sign and like-sign muon pairs with the J/ identification cuts 
described in section 3.4. 

Assuming the same pt distribution for unlike-sign and like-sign muon 
pairs from background, 



T I uncorr 



+- 



K 



++ 



{p 



N. 



•T)+- 



++ 



N. 



{p 



+- 



1.60 ± 0.17 GeV/c 



is obtained where a suffix H — (++) represents unlike (like)-sign pairs. 
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count 


(pt) 


Mpt) 


unlike-sign pairs 


41 


1.59 GeV/c 


0.15 GeV/c 


like-sign pairs 


5 


1.51 GeV/c 


0.31 GeV/c 



Table 8: Average pt values and their statistical errors for unlike-sign and 
like-sign muon pairs with no corrections. 

3.7.2 Efficiency corrections 

The Pt dependence of the J/ijj detection efficiency has been obtained 
which is relatively small (shown in Fig. 45). The correction factor C^ff = 
1.06 is obtained with a 3% statistical error because of slightly higher efficiency 
for low Pt J/ip^s. 

A systematic error of Ce// has been estimated to be ±4%, by increasing 
efficiency for lower pt {pt < 2 GeV/c) J/i/j^s while decreasing for higher pt 
(Pt > 4 GcV/c) by statistical errors of the simulation and vice versa. 

It is confirmed that a correction of (pt) due to rapidity dependence of 
efficiency is small (1%) since the correlation between pt and rapidity is small. 

3.7.3 Momentum offset and resolution correction 

The change in (pt) from finite momentum resolution and momentum shift 
has been estimated using a simulation, based on the observed center position 
of the J/tjj mass, 3156 ± 74 MeV/c^, and its resolution, 257 ± 75 MeV/c^ 
which are shown in Fig. 23. 

Assuming the shift of momentum is a constant (momentum-independent) 
and the effect of momentum resolution on (pt) is small, momentum shift is 
naively expected to be (3156+74-3097)/3097 ~+4% at the worst case (using 
Mj/^ = 3097 MeV/c^ [62]). To estimate the entangled effect of both the 
momentum shift and resolution on (pt), following two models are assumed 
for the shift 5p/p and resolution Ap/p: (I) Sp/p — a, Ap/p = b and (II) 
Sp/p = a'p, Ap/p = b'p, where a, a', b and b' are parameters which have been 
determined to reproduce the real data. For example, mis-calibrated magnetic 
field can cause (I), whereas mis-alignment of the MuTr chambers causes (II). 
The shifts of (pr) are estimated using PYTHIA events whose Pt shape is con- 
sistent with the real data as shown in the next section. The correction factor 

Cmom-offset IS Calculated as Cmom-offset = {Pt) generated/ {Pt) smeared whcrC 

(Pt) generated IS the average Pt of unpolarized J/ip^s generated by PYTHIA 
and {pt) smeared IS the average pt of J/V^'s reconstructed from the smeared- 
momentum muons with the formula (I) or (II) above. In the worst case, where 
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parameters are fixed to reproduce the 3230-MeV/c^ center position and 250- 

MeV/c^ resolution of J/ip, Cmom-offset = 0.95 is obtained for (I) and 0.94 for 
(II) with 2% statistical errors of the simulation both of which are consistent 
with the naive expectation. Without mass shift, where parameters arc tuned 
to the 3097-MeV/c^ center position and 250-MeV/c^ resolution, (pr) docs 
not alter over the statistical errors of the simulation (2%) for both (I) and 
(II) cases. Consequently, the correction factor Cmom-offset — 0.97 ± 0.05 is 
obtained. 



4 Results and Discussions 

4.1 Differential and total cross sections 

Results of rapidity-differential and double-differential cross sections for 
J / ip production measured in the South Muon Arm are summarized in Table 9 
and 10 respectively with each factor for them shown in the equation (1). 
Systematic errors of each factor are summarized in Table 11. 



y, dy, Pt, dpr 


= Nj/. 


^tot 

{%) 


C 

(nb-i) 


Brda/d,y (nb) 


1.70, 1.0, all, - 


41 - 5 = 36 


1.21 


81 


37.4 ± 7.1 (stat.) ± 10.1 (syst.) 


1.45, 0.5, all, - 


29 - 3 = 26 


1.31 


81 


49.1 ± 10.8 ± 13.3 


1.95, 0.5, all, - 


12 - 2 = 10 


1.08 


81 


22.9 ± 8.5 ± 6.2 



Table 9: Results of rapidity-differential cross-sections and their input values. 



y, dy, pt, dpr 


7V+_ - 

= Nj/^, 


^tot 

(%) 


C 

(nb-i) 


Br ( 1 / 27ipT)da / dpxdy 
[nb/(GcV/c)2] 


1.70, 1.0, 0.5, 1.0 


12 - 2 = 10 


1.40 


81 


2.81 ± 1.04 (stat.) ± 0.76 (syst.) 


1.70, 1.0, 1.5, 1.0 


17-2 = 15 


1.21 


81 


1.63 ± 0.47 ± 0.44 


1.70, 1.0, 2.5, 1.0 


9-1 = 8 


1.14 


81 


0.55 ± 0.22 ± 0.15 


1.70, 1.0, 4.0, 2.0 


3-0 = 3 


0.92 


81 


0.08 ± 0.05 ± 0.02 



Table 10: Results of double-differential cross-sections and their input values. 
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Vaiial)les 


Syst. errors (%) 




10 




A {Jj^ pol.) 


10 


J/i' 


^MuID 


11 


J/^ 
'^MuTr 


10 




Jl^ 
^BBC 


10 




Nj/^ 


negligible 


C 


^inela 


3 




^BBC 


15 


Total 


27 



Table 11: Summary of the systematic errors of each factor for the differential 
cross sections. 

4.1.1 Transverse momentum distribution 

Transverse momentum (pt) distribution for J/i/j production cross-sections 
at low pt {pt < Mj/^), which dominates the sample obtained, is sensitive to 
the intrinsic parton transverse-momentum (kx) as well as to initial momen- 
tum distribution of partons (mainly gluons) , but not to the production mech- 
anism much. With an increased yield especially at higher Pt {pt ^ Mj/^), 
we can discuss discrimination between the production models based upon 
a Pt shape. Here our data is compared with a phenomenological function 
which has been used to fit lower-energy results as well as a PYTHIA predic- 
tion based on the color-singlet model. 

Figure 51 shows invariant (double-differential) cross sections for J/i(j pro- 
duction as a function of pt- Error bars include statistical errors only. Lines 
show the function form 



da 

oc 



l+(^ 



with different values of a parameter B, which is related to the average pt 
{(Pt)) with (pt) = (1057r/768)i? by a numerical integration. All the lines are 
normalized to the data. This function form is consistent with our data with 
(pt> in the range 1.5 < {pr) < 2.2 GeV/c. It also fits the FNAL-E672/E706 
data [39], shown in Fig. 52, and the FNAL-E789 data [37] at y/s = 38.8 
GeV well. The histograms in those figures show the PYTHIA distributions 
which are consistent with both the real data and the fits with appropriate 
parameters.^ Contribution of 6-quark decays is neglected, since it is estimated 

''The low-pT cut off for partonic subprocesses is deactivated (default = 1 GeV/c). 
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to be small (1%) at ^/s — 200 GeV from other experiments and simulation, 
which will be described in Appendix A. The transverse momentum pt of a 
J/ip is sensitive to the rms of kr, (kr), which is set to 1.5 GeV/c at ^/i = 
200 GeV and 0.82 GeV/c at ^/s = 38.8 GeV to reproduce the real data as 
shown in Figs. 51 and 52. The result of (pr) will be discussed in the next 
section. 




Figure 51: Invariant cross sections for the J/ijj production as a function 
of Pt- The error bars include statistical errors only. The lines shows the 
function form da/dpr oc [1 + with different values of the average 

Pt, (pt) = (1057r/768) B GeV/c. All the lines arc normalized to the data. 
The histogram is a PYTHIA prediction with (/cr) = 1.5 GeV/c. 



4.1.2 Rapidity distribution and gluon distribution 

Since the gluon-fusion subprocess dominates at high energies, rapidity 
distribution for J/ip production cross-sections, da/dy, is sensitive to gluon 
distribution g{x, Q) in the proton, where a; is a momentum fraction of a gluon 
(Bjorken-x) and Q is a scale for the partonic interaction, which is in the order 
of Mc (charm quark mass) for the case of J/ip production. Sensitive x-region 
probed by our measurements is discussed, followed by the comparison of our 
data with some typical gluon distributions currently available. 

For low-pT J/ip production (p^ < Mj/^), da/ dy is roughly proportional to 
the product of gluon distribution function g{xi)g{x2) where Xi(2) is Bjorken-x 
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Figure 52: Differential cross sections for J/ip production as a function of 
at ^/s = 38.8 GeV. Data points are taken from the FNAL-E672/706 data 
in p+Be collisions [39]. The line shows the fit of the function da/dp\ oc 
[1 + (pt/-B)^]~*^ to the data. The histogram shows a PYTHIA prediction 
with {kr) = 0.82 GeV/c. 



of each gluon. They are related with y (rapidity) of a J/ if) by 

X\ — -\/^exp(y) and X2 — y/r exp{—y) 
where r = [Mji^j ^fsf' . The relation of X\ and for J jil) production events, 

X\Xi — T oc - (7) 
s 

is shown in Fig. 53 using PYTHIA and the GRV94LO PDFs with a scale 
Q — Mc — 1.45 GeV. The clear correlation XiX2 = r ~ (0.015)^ is seen with 
a band which represents finite px of J/t/j^s. The sensitive regions of x are 
0.01 < a; < 0.05 for J/^'s measured in the Central Arms and 2 x 10~^ < 
X < 0.01 and 0.05 < a; < 0.2 measured in the Muon Arms, which are also 
shown in Fig. 54 together with some gluon distributions of some typical 
PDF sets with a scale Q — Mc — 1.45 GeV. Distributions of Feynman-x 
{xf — xi —X2) of lower-energy experiments, obtained as xp — '^Pz/Vs where 
Pz is the z-component of the momentum of a J/ip, are well reproduced by 
these gluon distributions. For example, the xp distribution of the FNAL- 
E705 experiment at = 23.8 GeV [32] is weU reproduced by the GRV94LO 
PDFs as shown in Fig. 55. This agreement confirms that J/ip production 
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is dominated by the fusion of gluons, whose distribution is reasonable in 
the corresponding x-range, which is around a; ~ 0.1 covered also by the 
PHENIX Muon Arms as shown in Fig. 54. Another J/ip event generator, 
the RV generator [66], gives consistent xp distribution with PYTHIA, which 
confirms that xp distribution (or rapidity distribution at low pt) depends 
on gluon density but not on the event generator with a specific production 
mechanism. 




Figure 53: The correlation of xi and X2 of gluons probed by J ftp measure- 
ments with the Central Arms and Muon Arms. The line shows the low-pT 
limit, where xiX2 = r ~ (0.015)^. 



Using these PDFs, our rapidity distribution is well reproduced as shown 
in Fig. 56 including the result of the J/ip e+e~ channel measured in the 
Central Arms {da/dy \ y=o = 52 ± 13 ± 18 nb) [52]. Two curves in the figure 
show PYTHIA predictions with the GRV94LO and CTEQ4M PDF sets fitted 
to the data, both of which are consistent with the results. Their difference 
is small compared to the errors of our measurement. Other PDF sets are 
also tried and fluctuation is found to be small (3%) in the determination of 
CTj/^ described next. This consistency confirms reasonable gluon density in 
the x-region 2 x 10"^ < x < 0.2. 
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Figure 54: The gluon density xg{x) of three typical PDF sets with a scale 
Q — Mc — 1.45 GeV, together with sensitive x regions probed by J/t/j mea- 
surements with the PHENIX and FNAL-E705 experiments. 
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Figure 55: Feynman-x (xp) differential cross sections for J/ip production 
obtained by the FNAL-E705 experiment in p+Li collisions at ^/s — 23.8 
GeV [32]. The hne shows a PYTHIA prediction with GRV94LO PDFs. 




Rapidity 



Figure 56: Rapidity differential cross sections for J/ip production. Both 
statistical and systematic errors are shown, which are in bars and in brack- 
ets respectively. Two lines show PYTHIA predictions with different parton 
distribution functions (GRV94LO and CTEQ4M) fitted to the data. 
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4.1.3 Total cross section and its compcirison with lower-energy 
results 

Assuming a PYTHIA rapidity distribution which reproduces our results 
well and the common branching fraction Bu for the electron decay channel 
(Bee) and the muon decay channel (Bf^fj), 

Biiaj/^ = 226 ± 36 (stat.) ± 79 (syst.) nb 

is obtained. Systematic error is obtained as the maximum plausible system- 
atic spreads including both the electron and muon measurements. Using the 
average value of B^e and B^^ ^ for Bu, total production cross-section 

c'j/V' — 3-8 ± 0.6 (stat.) ±1.3 (syst.) //b 

is extracted. 

Figure 57 shows ^/s dependence of ctj/^ together with predictions using 
gluon densities of three different PDF sets with a scale of Q = M^. Lower- 
energy experimental results are summarized in Table 12. Our result of the to- 
tal cross section is at the highest energy since higher energy {^/s > 630 GeV) 
experimental data on J/ip^s suffer from a Pt cut {px > 4 or 5 GcV/c) as de- 
scribed in section 1. The increase in aj/^ at higher y/s can be explained by 
the increase in gluon density in lower x-region according to the relation (7). 
For the absolute normalization for aj/^, color-octet matrix elements are used, 
whose uncertainties will be described later. Charm quark mass is adjusted 
to Mc = 1.45 GeV/c^ or = 1.50 GeV/c^ depending on PDF sets to re- 
produce the normalization best. Parton distributions are sensitive to a scale 
Q. Figure 58 shows dependence of cxj/^ with different values of Q using 
the GRV94LO PDFs and = 1.45 GeV/c^. For higher Q, distribution is 
steeper since gluon density increases more in lower-x than in higher-x as Q 
increases. The experimental data prefer Q = Mc for this case. In spite of 
uncertainty from the choice of a PDF set and a scale, ^/s dependence of cxj/^ 
including our new result is consistent with gluon distribution in the proton. 

The absolute normalization for aj/^ depends on the production models. 

A consequence of the GEM is the constant ratio, that is ^/s and process 
independent, of aj/.^ to ci^)/) where o"/)/) is the total production cross section 
for D mesons. It is determined from the photo-production data to be about 
0.06 [14]. To normalize a^jj in hadro-production, the phenomenological K- 
f actors are used to account for higher order corrections. As a result, crj/^ of 
hadro-production data are well reproduced including our measurement. 

*^5ee = 0.0593 ± 0.0010 and = 0.0588 ± 0.0010 [62] 
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Target 


(GeV) 


aj/^ (nb) 


(Pt) (GeV/c) 


Reference 


P 


6.7 


0.62 ±0.18 


— 


[26] 


Be 


7.3 


9+2.0 
^-1.0 




[1] 


P 


8.6 


2.4 ± 1.2 




[27] 


Be 


11.5 


22 ± 6 




[28] 


Be 


16.8 


138 ± 46 


0.96 ± 0.12 


[29] 


P 


16.8 


94 ± 20 




[30] 


P 


19.4 


122 ± 22 


0.95 ± 0.02 


[30] 


C 


20.5 


190 ± 26 


0.98 ± 0.04 


[31] 


Li 


23.8 


162 ± 22 




[32] 


P 


24.3 


144 ± 19 


1.1 ± 0.2 


[34] 


Be 


27.4 


220 ± 54 




[35] 


P 


30 




1.14 ± 0.12 


[42] 


Be 


31.5 


322 ± 70 


1.08 ± 0.05 


[36] 


Be 


31.5 




1.15 ± 0.02 


[39] 


Be 


38.8 




1.22 ± 0.01 


[39] 


P 


52 


700 ± 320 


1.0 ± 0.2 


[40] 


P 


53 




1.39 ± 0.05 


[42] 


P 


63 




1.29 ± 0.05 


[12] 



Table 12: The J/ip total cross sections aj/^ and the average values of pt 
{{Pt}) in proton-induced interactions at lower energies. Target mass depen- 
dence is accounted for assuming A'^-^ [3] where A is the atomic number of 
the target. 
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Figure 57: Center of mass energy dependence of the production cross section 
for J 1"^ in nucleon-nucleon collisions. Lower energy results are summarized 
in Table 12. Error bars include both statistical and systematic errors which 
are added in quadrature. The curves are for QCD predictions with different 
PDF sets with a scale Q — Mc- 




Figure 58: Center of mass energy dependence of the production cross sec- 
tion for J/ip in nucleon-nucleon collisions with different PDF scales [Q). 
GRV94LO PDFs and charm quark mass Mc = 1.45 GeV/c^ are used. 
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The CSM (color-singlet model) has been known to fail to describe the 
total J/ip cross sections at lower energies by a factor of about 20 [22]. Despite 
large theoretical uncertainties from the choice of scale parameters and charm 
quark mass, this large discrepancy can be hardly explained. 

Color-octet contributions have been considered [22-24] in 2 — > 1 subpro- 
cesses (dominantly g + g J/i^) which is the largest contribution to the low 
Pt (or pr-integrated) yield. Cross section for the direct J/ip production with 
gluon-gluon fusion can be written as 



X 



where is the strong coupling as a function of a rcnormalization scale (^u). 
The color-octet matrix elements {O'^^^ i^^^^Lj)) , each of which represents 
probability for a color-octet cc pair in a state "^^^^Lj to form a, J/ if) meson, 
have been extracted from experimental data [23] . There are also decays from 
the color-octet production of Xci but their contribution is neglected, since it 
is estimated to be small compared to the color-singlet production of Xc [23]. 
Including this color-octet contribution to the direct J/tfj production, total 
cross section agrees with the data as shown in Figs. 57 and 58 using reasonable 
charm quark mass and a renormalization scale The uncertainty of uj/^ 
from Afe, which is not constrained well yet [62], has been studied by changing 
it by about 10% and variation of ctj/^ is found to be about a factor of two. 
The effect of the uncertainty of ii is also studied by making it double or half 
and (Tj/^ changes by a factor of two to three. The normalization for aj/^ with 
the COM is consistent with the experimental data within large theoretical 
uncertainties inherent in perturbative QCD. 



4.2 Average px 

The average value of px of {pr) is sensitive to the intrinsic par- 

ton transverse-momentum (/c^) as well as initial momentum distribution of 
partons (mainly gluons), but not to the production mechanism much. It 
is calculated as (pr) = CeffCmom-offset{PT) uncorrected where each factor has 
been obtained in section 3.7. Table 13 summarizes their statistical and sys- 
tematic errors. As a result. 



(Pt) = 1.66 ± 0.18 (stat.) ± 0.12 (syst.) GeV/c 
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is obtained with cuts < pt < 5 GeV/c and 1.2 < y < 2.2. Correction 
factors on {pt) to extrapolate it to the entire kincmatical region are estimated 
in the next, followed by the comparison of {pr) with lower-energy results. 



Statistical 


± 11% 


Efficiency corrections 


± 5% 


Momentum offset 


± 5% 


Syst. total 


±7% 



Table 13: Summary of statistical and systematic errors of (pr)- 



High-pT correction 

Since the statistical hmit of the data is < 5 GeV/c, the effect on 
{Pt) from missing high-p^ events is estimated assuming the function form (6) 
which fits the FNAL-E672/E706 data [39] (Fig. 52) the FNAL-E789 data [37] 
well. As a result, (pr (0 < Pt < 5 GeV/c)) = 1.82 GeV/c and {pr {0 < Pt < 
ooGeV/c)) = 1.88 GeV/c are obtained. Therefore a correction factor from 
high-pT events is expected to be small (3%). 

Other two empirical functions forms which have been used for other ex- 
periments, 

^ (X exp(-apT) (8) 

and 

cx exp(— 6p^) (9) 



2 



dpi: 

arc tried where a and b are fit parameters. Their fits to OTir data is shown 
in Fig. 59. For (8), pr slope is slightly harder than other functions but still 
consistent with our results. The correction factor of 11% is obtained. For 
(9), much smaller a correction factor, 0.04% is obtained. The reduced chi- 
squares for the fit of (8) to the FNAL-E789 data (xV^DF = 6.6) is much 
worse than others (xV^DF = 0.5 for (6) and 1.3 for (9)). 

Consequently, high-p^^ correction (p^ > 5 GeV/c) is expected to be small 
(< 3%) assuming the function form for the px spectrum which fits our results 
as well as the FNAL-E789 and FNAL-E672/706 results at = 38.8 GeV 
with good chi-squares. 
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Figure 59: Invariant cross sections for the J/ijj production as a function of 
Pt together with three different function forms. 



Rapidity correction 

Figure 60 shows rapidity dependence of (pt) using PYTHIA with the 
GRV94LO PDFs and other parameters which reproduce px and rapidity 
distributions of both lower energy experiments and our result as already 
shown. Using this dependence {pt) i.i<y<i.2 / {pt) -oo<y<oo is found to be 0.99 
which is close to unity. 

Comparison with lower energy results 

At higher energies, more energetic partons participate in interactions, 
thus causing a slight increase in (p^) of J/ip^s. Figure 61 shows \/i de- 
pendence of (pt) including the results of PHENIX and the lower energy 
experiments summarized in Table 12. Neither high-p^ correction nor rapid- 
ity correction is applied to the PHENIX data, which are estimated to be 
small, as well as the other experiments. Our result, at the highest energy, 
is slightly higher than the lower-energy results, which is consistent with the 
expectation. The experimental data are well fitted with a function p+g In ^/s 
with a good reduced chi-square (x^/NDF ~ 1.1) as shown in the same figure, 
where p and q are fit parameters. On the other hand, another empirical form 
p' + which fits the lower energy results well [36], does not fit including 
the PHENIX result, indicating it is applicable only in the limited ^/s range. 
Another curve in the figure shows the PYTHIA prediction with a fixed {kr) 
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Figure 60: Rapidity dependence of {pt) in y/s = 200 GeV p+p collisions 
obtained with PYTHIA and the GRV94LO PDFs. The hatched area shows 
the South Muon Arm acceptance. 



= 0.82 GeV/c, which reproduces the pt spectrum of the FNAL-E672/E706 
data at ^/s = 38.8 GeV as shown in Fig. 52. The PYTHIA prediction is 
consistent with our data, but slightly higher than the experimental results 
at lower energy (y^ < 20 GeV), indicating (kr) increasing with y/s [67,68]. 

Such ^/s dependence has also been observed in (pt) of muon pairs from 
the Drell-Yan process, which suffers less theoretical uncertainties such as 
next-to-leading-order corrections thus has direct relation to (kx) [68]. Its 
magnitude is also consistent with that of J/ip at lower energies [37,39,69, 
70]. This similarity supports that (pt) of J/ip can be understood by a 
perturbative-QCD based model (PYTHIA) with reasonable (/ct) values. 



4.3 Future measurements 

Our result on the total cross section for the J/ip production in = 
200 GeV p+p collisions is consistent with both the color-evaporation model 
and color-octet model predictions. In this section, we will discuss future 
measurements expected at RHIC which will provide further critical tests for 
these models. 
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Figure 61: Center-of-mass energy dependence of average pr of J/ip^s includ- 
ing the results of PHENIX and lower energy experiments. The error bars 
include both statistical and systematic errors which are added in quadra- 
ture. The solid line shows a fit with a function p + qln ^/s where p and q 
are fit parameters. The dotted line shows the PYTHIA prediction using the 
GRV94LO PDFs and parton-transverse-motion (/c^) = 0.82 GeV/c which is 
tuned to reproduce the px spectrum of the FNAL-E672/E706 data at ^/s — 
38.8 GeV. 
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Polcirization (Spin Alignment) of J/ip 

One critical test for these models is measurement of polarization A (or 
spin-alignment) of the J/t/j. The GEM predicts zero-polarization due to 
multiple gluon emissions while the CSM and COM predict sizable polariza- 
tion. All the lower energy experiments show zero-polarization (|A| < 0.1) [43] 
which is consistent with the GEM prediction except at large xp (A — 1 
with Xp 1). The CDF experiment measures polarization of J/ijj^s with 
Pt > 4 GeV/c [47]. Medium-pr J/^/^'s (4 < pr < 12 GeV/c) are transversely 
polarized (A > 0) and consistent with the COM prediction. However, at 
high pt {pt > 12 GcV/c), polarization decreases with pt unlike the COM 
prediction, causing a sizable discrepancy at the highest pr bin (15 < pr < 
20 GeV/c). 

Our measurement of A at y/s — 200 GeV will serve as another vital input 
for the production models. Since A is obtained by fitting cos 9* distribution 
to the form 

"^"^ ocl + A(cosr)^ 



d cos 9* 

sensitivity to higher | cos^*| is important for the measurement of A, where 9* 
represents the decay angle of a muon in the J/ip rest frame with respect to the 
momentum direction of the J/ijj in the laboratory frame (Gottfried- Jackson 
frame) . 

Figure 62 shows the acceptance of the Muon Arms for J/ip^s with px > 
2 GeV/ c, as a function of | cos 9* \ . Our acceptance is sensitive up to | cos 9* \ ~ 
0.8. Acceptance for larger |cos^*| increases with pp of J/ip^s, because the 
backward-going muon in the J/ijj rest frame is boosted toward the detector 
in the laboratory frame. 

Polarized cross section for the J/ip production 

In Run-2, RHIC has successfully accelerated transversely polarized pro- 
tons up to 100 GeV/c momentum with about 20% polarization and coUided 
them at ^/s — 200 GeV. Details of spin physics at RHIC will be found in [8]. 

In Run-3 longitudinally polarized p+p collisions with about 50% polar- 
ization are expected at RHIC. We plan to measure double-longitudinal spin 
asymmetries for J/ijj production {Aj^^^^ defined as 

^LL - (++) (+_), 

where Cj^^'* (cj/^ ■*) stands for production cross-section when incident protons 
have the same (opposite) signs of helicities. 
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Figure 62: The Muon Arm acceptance for J/ip^s with px > 2 GeV/c as a 
function of I cos 6* I . 



It is sensitive to polarized gluon density {/S.g{x) — g'^{x) — g~{x)) in the 
proton where g^^~\x) denotes density of gluons which are polarized in the 
same (opposite) direction as the proton helicity and can be written as 



where Xi and X2 are momentum fractions of each gluon. The partonic asym- 
metry a^^^^'^^^^ for the subprocess g -\- g ^ J/i/jX is calculated by each 
model in [71]: 

oll ~ +1 with the GEM, 
clll ~ —1 with the CSM and 

5M|e + l5B„{OrePo)) - iB,{0i"fP2)) 



where 



e = (03^/^(%))-^(03^/n^Po))and 
Bj^Br{xcJ^J/i^l)- 
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Calculation with the COM is for the 2^1 subprocesses, which dominate the 
low-pT yield. The mere determination of the sign of A"^^^ can choose either 
the CSM or CEM. For the COM, G is currently unknown and expected to 
be extracted from Aj^^ of photo-production data to be measured in polar- 
ized proton and electron collisions at the Electron Ion Collider (EIC) [72]. 
Using B obtained from the unpolarized photo-production data [73], branch- 
ing fractions [62] and the color-singlet matrix elements {O^"'' i^Pj)) from the 
data of hadronic decays of Xc states [74], with the color-octet model is 
expected to vary from —0.3 to 0.7 assuming two extreme cases, {O'l^'^ {}Sq)) 
= and {Og^^ {^Pq)) = 0. For the color-octet 2 — 2 subprocess, A'l^^ has 
also been calculated in [75] assuming a function form of Ag{x) found in [76], 
which is shown in Fig. 63 as a function oipr of J/tp^s together with statistical 
sensitivities of PHENIX expected with one RHIC-year luminosity.^ Measure- 
ment of Aj^^^ is motivated not only by the direct determination of polarized 
gluon distribution but also by another test for the production models for 
charmonia. 

0.14 
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0.1 
0.08 
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0.04 
0.02 


-0.02 
-0.04 

2 4 6 8 10 12 14 

p^(GeV) 

Figure 63: Double- longitudinal spin asymmetries for the J/ip production 
(Ajl^) as a function of p^^^ with COM expectations together with statistical 
sensitivities of PHENIX with a RHIC one-year full luminosity. 




^In the enhanced mode. Also in this paper, 100% detection efficiency for J/tp's are 
assumed. 
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Table 14 shows the summary of model sensitivities of each observable 
measured at RHIC. Our measurement of aj/^ prefers the GEM and COM 
rather than the CSM. Measurements of rapidity and low-pr distributions 
{{Pt)) do not have sensitivities to the models. Future measurements of the 
polarization of J/ijj and A'^^^ will provide discrimination of the GEM from 
the GOM. 





GSM 


GEM 


GOM 


comments 




X 





done 


rapidity distribution 


no sensitivity 


done 


(Pt) 


no sensitivity 


done 


Polarization 


> 





> 


needs high-p^ events 




< 


> 


vary 


needs ^g{x) 



Table 14: Summary of model sensitivities of each observables regarding J/ if: 
measured at RHIG. Our measurement of ctj/^ prefers the GEM and GOM 
rather than the GSM. Future measurements of the polarization of J/V' and 
^LL ^^^^ provide discrimination of the GEM from the GOM. Expected ranges 
of values for those measurements are shown with each model. 



5 Conclusion 

Measurements of J/'4> mesons in p+p collisions at RHIG provide useful 
information on both the perturbative and non-perturbative aspects of QGD. 
They also play crucial roles in both QGP physics and spin physics expected 
at RHIG. 

With the South Muon Arm in the PHENIX detector covering 1.2 < y < 
2.2, J/il^ particles have been successfully triggered by the simple trigger sys- 
tem and clearly identified with a small background via /x^/x" decays in the 
first p+p collisions at y/s = 200 GeV. Total and differential cross sections 
for inclusive J/ip production have been obtained. 

Transverse momentum distribution is consistent with the various function 
forms and PYTHIA which fit the other experimental data. Average trans- 
verse momentum {px) — 1-66 ± 0.18 (stat.) ± 0.12 (syst.) GeV/c is slightly 
higher than the lower energy results and consistent with the PYTHIA 's ex- 
pectation with a reasonable value of (/ct)- 

Rapidity-differential cross section for inclusive J/ijj production times the 
branching fraction to a pair B^^da / dy\y=i,-j = 37 ± 7 (stat.) ± 11 
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(syst.) nb has been obtained. Including the measurement at mid rapidity 
(|y| < 0.35) using the electron decay channel, rapidity distribution is found 
to be consistent with the gluon distribution function of some typical PDF 
sets in the Bjorken-x range 2 x 10~^ < a; < 0.2. 

The total cross section cxj/^ = 3.8 ± 0.6 (stat.) ± 1.3 (syst.) nh has 
been extracted from the fit of the rapidity distribution, which is related 
with the lower energy results by a gluon distribution function with a rea- 
sonable scale. Proper normalization for crj/^ is reproduced well by both the 
color-evaporation model and the color-octet model using reasonable gluon 
densities, scales and charm-quark mass. 

Our results confirm the validity of the production mechanism of J/ip 
based on perturbative QCD. The further studies of the production mecha- 
nism, via measurements of such as the spin alignment of J/ip and A'l^^ , are 
clearly needed especially to understand the non-perturbative aspects involved 
in the production of J/ip. The results have also established the baseline nec- 
essary for the QGP search in Au-|-Au collisions as well as gluon polarization 
measurements in polarized p-|-p collisions at RHIC. 

A Estimation of various contributions to the 
inclusive J/ip yield 

J/ijj mesons can be produced either directly or via decays from higher 
state charmonia (xc and ip') as well as from 6-quarks. Our measurement is 
inclusive and it is not possible to separate each contribution from others. 
We will estimate fractions of each contribution taking into account results of 
both higher and lower energy experiments. 

6-queirks 

Contribution from 6-quark decays is energy-dependent. Production cross 
section for a bb pair times the branching fraction for a b (or 6)-quark decaying 
into a J/ip meson in p (or p)+N collisions o'(p(p)N bbX)Br{bb J/ipX) 
has been measured in p+Au collisions at y/s = 38.8 GcV by the FNAL-E789 
experiment [38] and the result is 148 ± 34 (stat.) ± 28 (syst.) pb/nucleon. 
It has been also measured in p-|-p collisions at ^/s — 630 GeV by the GERN- 
UAl experiment [44] and 32.7 ± 2.6 (stat.) ± 13.4 (syst.) nb has been 
obtained with p^^^ > 5 GeV/c and \y''^'^\ <2.0 cuts. The ratio of those 
results is consistent with a perturbative QCD prediction using PYTHIA. It 
is extrapolated to = 200 GeV and a{p+p bbX)Br{bb J /i^X) = 40 
nb is obtained which is about 1% of our inclusive result. 
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Production cross sections for ip', or ip (2S) mesons, have been also mea- 
sured by the FNAL-E789 [38] experiment. The relative production cross 
section times the branching fraction to a muon pair for a ijj' to that for 
a J/ip, R = Briijy i2+i2-)a^,/Br{J/ip ii+ii-)aji^ = 0.018 ± 0.001 
(stat.) ± 0.002 (syst.) has been obtained. This ip^ fraction can con- 
tribute to our J/ip measurement because of the poor mass-resolution (200 
to 300 MeV/c^) currently observed. Using the branching fractions Br{'4>' — > 
= (7.0 ± 0.9) X 10-=^ and Br{J/ip //+//-) = 5.88 ± 0.10 % [62], 
o^ijuji^ is determined to be 0.15 ± 0.03. Contribution of ip' J/ip fj,'^fj,~ 
to inclusive J/ip yield is therefore expected to 15% x 0.55 {Br{ilj' J/ipX)) 
~ 8%. The color evaporation model expects energy and process independence 
of the fractions of ip' and Xc to the inclusive J/ip yield. Slightly higher or 
consistent R values have been obtained by the UAl (R — 0.029 ± 0.010 
(stat.) ± 0.007 (syst.) with pr > 5 GeV/c and \y\ < 2.0 cuts) [44] and CDF 
{R = 0.033 ± 0.002 (stat.) toms (syst-) nb with pr > 5 GeV/c and \r]\ < 
0.6 cuts) [45] experiments. 

Xc 

The fraction of the inclusive J/ ip yield from radiative x decays has been 
determined to be (30 ± 4) % by the FNAL-E705 experiment in p-|-Li coUisions 
at = 23.8 GeV [33]. A consistent value, (29.7 ± 1.7 (stat.) ± 5.7 (syst.))% 
was obtained in p-|-p collisions at ^/s = 1.8 TeV by the CDF experiment at 
Fermilab [46]. 

In summary, the fractions of 6-quark, ■0' and x contributions to inclusive 
J/ijj yield at RHIC energy arc expected to be 1%, 10-15% and 30% respec- 
tively. The remains are ascribed to the direct production. The fractions 
of ip' and Xc contributions do not agree with the color-singlet model predic- 
tion. Including color-octet contributions, they are successfully explained [23]. 
The color-evaporation model can not predict their fractions. However, small 
^/s dependence of them in nucleon-nucleon collisions is consistent with its 
assumption that they are process- independent. 
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B Attached figures' captions 

• figure 1: The RHIC accelerator complex, 
figure 2: The PHENIX experiment overview. 

figure 3: Magnetic field lines inside the Central and Muon Magnets. 

figure 4: A photograph of a Beam-Beam Counter (BBC). Dimensions 
are 29 cm in outer-diameter and 25 cm in length. 

figure 5: A schematic view of a Muon Arm. It consists of a conical 
magnet (Muon Magnet) with tracking chambers inside at three stations 
(Muon Tracker, or MuTr) and an array of chamber planes interleaved 
with absorbers (Muon Identifier, or MuID). 

figure 6: Cross section of a MuTr station. One MuTr station consists of 
two or three gaps, each of which has one anode-wire and two cathode- 
strip planes. 
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figure 7: Octant and half-octant structures of a MuTr gap. The beam 
hne passes at the center of this figure. The bold lines show the bound- 
aries for octants and this lines for half-octants. Octant and half-octant 
numbers are also shown. For the South Arm, the left-hand side in the 
figure (octant-5 side) corresponds to the cast direction. 

figure 8: A photograph of a MuTr station-2 octant. A cathode plane 
is made of etched 25-/xm copper coated mylar foils. 

figure 9: Directions of cathode strips and anode wires in station- 1, 
octant-2, half-octant-0 planes. Bold lines show cathode strips in non- 
stereo-angle planes. Dotted lines show cathode strips in stereo-angle 
planes with angles shown in Table 4. Angles in the figure are not exact. 

figure 10: The optical alignment system in which light from an optical 
fiber is projected from station 1 through station 2 to a CCD mounted 
on station 3. The relative chamber positions arc monitored to ±25 /im. 

figure 11: A schematic diagram for the MuTr Front-End Electronics 
(FEE). Raw chamber signals are continuously amplified with CPAs 
(Charge Pre- Amps) and stored in AMUs (Analog Memory Units). Upon 
receipt of a level- 1 trigger bit from a GTM (Granule Timing Module), 
stored samples from all channels are digitized by ADCs (Analog to Dig- 
ital Converters) and the results are sent to a DCM (Data Collection 
Module. 

figure 12: Measurement of position resolution obtained in cosmic-ray 
tests of a station-2 octant. The composite chamber plus projection 
error (position resolution for a track obtained with the other five planes) 
was about 131 fim, consistent with the 100 /im specification for the 
chambers and readout alone. 

Figure 13: Setup for the MuID beam test carried out at KEK-PS. Four 
scintillation counters (STl to ST4) define the beam. Three gaseous 
Cerenkov counters (GCl to GC3) identify pions and muons. 

Figure 14: Pion mis-identification rate as a function of pion momentum 
for both the South and North Muon Arms obtained with the beam test 
at KEK. The solid (dotted) fines show simulation results including 
(excluding) tt^ — > /I'^u decays. 

Figure 15: Schematic diagram of the MuID NIM-logic trigger. Each 

"Quadrant" MLU takes 16 (4 for each gap) output signals from the 
MuID FEE (Front-End Electronics). The "Decision" MLU takes out- 
put signals of all the "shallow" and "deep" quadrant MLUs and gen- 
erates final trigger bits, which are IS (single shallow), 2S (double shal- 
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low), ID (single deep), IDIS (deep-shallow) and 2D (double deep), as 
PHENIX Global-Level- 1 input. 

figure 16: An accepted-pattern example of the MuID NIM-logic trigger. 
Capital letters represent positions of each segment. Gap 2 was not used 
because of the limitation of input channels of the trigger logic. 

Figure 17: The integrated luminosity RHIC has delivered to the PHENIX 
and STAR experiments as a function of day starting from December 
21, 2001. 

Figure 18; BBC efficiency for p-|-p ir^X events as a function oi pt 
of measured in the Central Arms (|?7| < 0.35) [64]. 
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